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EDITORIAL 


At the time of writing, the score in the satellite-launching contest is two-three. 
Another Sputnik is ready for the May Day celebrations, the U.S. Navy plans 
to launch seven full-size, fully-instrumented satellites this year, and more 
ambitious versions of the Explorer are in the offing. 


However, while to the man-in-the-street these American and Russian 
programmes must sometimes take on the aspect of a (very expensive) game, 
they are nevertheless important scientific experiments. It is interesting to 
note that each of the five satellites which have so far been placed in orbit has 
had its own individual characteristics as regards shape, orbit and instrumenta- 
tion. One cannot draw any conclusion as to the relative value of the two sets 
of data obtainable from a “‘grapefruit’’ which stays in orbit for five years, and 
from keeping a live dog in orbit for a few days (or a dead one for months). They 
are complementary, and both are valuable. 


Research in outer space has come to stay; one astronomer has even become 
worried that a commonwealth of satellites flitting round the globe will interfere 
with his observations. But there is no true commonwealth at the moment: only 
two nations are participating. This state of affairs cannot continue indefinitely. 
The Duke of Edinburgh spoke for us all when he said recently that he could 
not believe that the people of the commonwealth would be content to sit by 
and watch others explore the universe. Sooner or later other nations will 
have to start programmes of space research, if only to bolster up their prestige 
as possession of the hydrogen bomb spreads to wider and wider circles. 


Has not the time come for a really international agency to be set up under 
the aegis of the United Nations, so as to enable scientists of all races to partici- 
pate in this work? An immediate objection to such a proposal is, of course, 
the question of military security, but is that objection valid? Three years 
ago it was found possible to declassify a considerable volume of information on 
nuclear energy, and no political or military débacle resulted. It is doubtful 
whether one would follow a similar release of rocket engineering data or the 
setting up of a U.N. Commission for Space Research. 
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INTERPLANETARY ORBITS* 


By M. VERTREGT,t MEMBER 





SUMMARY 


The basic equations under simplified conditions for interplanetary flight are derived. 
For a voyage from planet to planet an unlimited number of orbits is possible. In 
order to give a clear survey of these possible orbits a diagram is developed from which the 
approximate energy-requirement, the duration, and other particulars of a voyage can be 
easily found. 


LIST OF SYMBOLS 
A the position of the Earth at the moment of departure of the spaceship. 
a the semi-major axis of the orbit of the spaceship. 
E  ameasure for the required energy of an orbit. 
€ the numerical eccentricity of the orbit. 
@ =the angle between the radius vector and the tangent of the orbit. 
p __ the constant of the Sun’s gravitational field (u = 1-329 x 10** cm. sec.~*). 
n the ratio r,/7,. 
P, the position of the planet at the moment of departure of the spaceship. 
P, the position of the planet at the moment of arrival. 
q the ratio a/r,. 
vy, the radius of the orbit of the Earth. 
vy, the radius of the orbit of the planet P. 
vy, the perihelion-radius of the orbit of the spaceship. 
Y439 the aphelion-radius of the orbit. 
S the Sun. 
T the duration of the voyage. 
7 an angle used in calculating the duration of the voyage. 


¢, the angular perihelion-distance of the Earth on the orbit of the spaceship at the 
moment of departure. 


¢, the angular perihelion-distance of the planet on the orbit of the spaceship at the 
moment of arrival. 


@ the angle between 4 and P,, or the total angular distance of the voyage. 
¥ 


the angle between A and P,, or the configuration of the two planets at the moment of 
departure. 


’, the linear velocity of the Earth in its orbit. 
the linear velocity of the planet in its orbit. 


U 

J 

V, the linear velocity of the spaceship in its orbit at position 4. 
V, the linear velocity of the spaceship in its orbit at position P,. 
V 


‘iq the component of the velocity of the spaceship caused by the engine at the moment of 
departure. 


V,, the component of the velocity of the spaceship caused by the engine at the moment of 
arrival. 
1. Introduction 
In a future time spaceships will ply between the planets as in the present 
time our seaships and airplanes ply between the continents of the Earth. No 


* Manuscript received 14 March, 1957; Addendum received 30 October, 1957. 
t Mozartlaan 73, The Hague, Netherlands. 
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one can tell if this will be realized in a far or in a near future, but if we consider 
the rapid development of flying in the course of the past half-century a modicum 
of optimism is certainly justified. 

But before we can speak of real spaceflight certain conditions must be 
fulfilled. In the first place, we must have at our disposal more powerful 
sources of energy than are available now in the form of chemical energy. 

When we consider that nowadays the simplest task in spaceflight, namely 
the launching of a small unmanned artificial satellite, taxes the resources of 
two great nations and requires a payload ratio of about 1,000, we cannot 
expect to be able to cope with the much more exacting task of establishing 
communication between the planets with this insufficient source of energy. 

Therefore we must assume that in the future means will be found to adapt 
the vastly more potent source of nuclear energy to the propulsion of spaceships, 
and that characteristic velocities of at least 50 km. per sec. will be attained 
without excessive payload ratios. 

Only if this condition is fulfilled will we possess a reserve of energy in our 
spaceship which will enable us, if necessary, to change our velocity or our 
course—in short, to navigate a ship in space. 

Though at the present time this may sound utopian, there is nothing 
intrinsically impossible in this condition and we may safely assume that it will 
be realized at some time in the future. 





2. Navigation in Space 

Navigation on Earth consists of setting a course between two fixed points 
and keeping that course. In space there are no fixed points except the Sun, 
which we consider as such. Both our starting point and our point of destina- 
tion describe orbits around the Sun, and in the case of a departure from and an 
arrival at an artificial satellite they also describe orbits around the planets. 
For these reasons navigation in space is vastly more complicated than naviga- 
tion on Earth and special techniques must be used. 

In a previous paper’, the author described a method for orientation in space; 
in this paper we will consider some of the problems connected with navigation 
in space. 

In order to simplify as much as possible the rather complicated circum- 
stances, we make the following assumptions: 

(1) We assume that the greatest part of the voyage is made in free flight, 
i.e., that the duration of the periods of acceleration and deceleration are 
insignificant compared with the total duration of the voyage. 

(2) We assume that the planets do not exercise an attraction on the space- 
ship, but we take into account only the attraction of the Sun. As a 
consequence of assumptions (1) and (2) we can consider the orbit of a 
spaceship to be a true Kepler orbit. 

(3) The planets describe circular orbits around the Sun. 

(4) The planes of the orbits of the planets are situated in the plane of the 
ecliptic. 

(5) Only interplanetary voyages are considered. 
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(6) Unless specified otherwise, we take the Earth as the planet of departure. 
The equations given in this article can be used for any other planet with 
only small alterations. 





Fic. 1. Motion in an elliptical orbit between the Earth (A) 
and an inner planet (P). 


Fig. 1 shows an orbit of a spaceship between planet A (the Earth) and an 
inner planet P around the sun S. 
The angles in this figure are: 
¢@, the angular perihelion-distance of the Earth on the orbit of the spaceship 
at the moment of departure. 
¢, the angular perihelion-distance of the planet on the orbit of the space- 
ship at the moment of arrival. 
@ the angle between A and P,, or the total angular distance of the voyage. 
WY the angle between A and P,, or the configuration of the two planets at 
the moment of departure. 
The distances are: 
y, the radius of the orbit of A (the Earth). 
r, the radius of the orbit of the planet P. 
r, the perihelion-radius of the orbit of the spaceship. 
1199 the aphelion-radius of the orbit of the spaceship. 
The elements of the orbit of the spaceship, which can be an ellipse or a 
hyperbola, are: 
a the semi-major axis, 


e the numerical eccentricity. 
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In order to simplify our equations we express all distances in the length of 
the semi-major axis of the orbit of the Earth, the Astronomical Unit. 
Therefore we call: 


'o 


-= a - i re as (1) 
" 
a 
- =¢ oii - “ ea es (2) 
e 


3. Possible Orbits 

Between the points A and P, (Fig. 1) an unlimited number of ellipses and 
hyperbolae with focus S can be drawn. All these conics are theoretically 
possible orbits for spaceships, but they all require different energies, they have 
different durations, and the configurations of the planets at the moment of 
departure are different. It is the aim of this paper to get a clear survey of 
the particulars of these orbits. 

Firstly we will investigate the general conditions for an orbit that can be 
used for interplanetary flight. 

From Fig. 1 it is obvious that the perihelion-radius 7, of an orbit to an 
inner planet must be smaller than or at the utmost equal to the radius of the 
orbit of the planet, and furthermore that the aphelion-radius 7,.5 of the orbit 
must be greater than or at least equal to the radius of the orbit of the Earth. 


Therefore we can put for inner planets: 


%, <1%e Yun > 1, 7 yf ae (3) 
We know that: 

% = a(l — €) Yigg = (1 + €) 
Thus: 

a(l—e)<% a(l+e«) >”, of = bf (4) 
Dividing by r, we get: 

gil—e) <n gi+e)>1 nes ms ee (5) 
From which: 

np} Boe € ee in iia sul (6) 

q q 


For a hyperbolic orbit we have only the limitation that for an inner planet the 
perihelion-radius 7, must be smaller than or at the utmost equal to the radius 
of the orbit of the planet. 


Therefore: 
ae - ee Me sa - (7) 
%=ale—1) .. im a bo 3% (8) 
or: 
ale—I1) <n, <> _ _ aoe - (9) 


Dividing by 7,: 
qie—1)<n... - * i .. (10) 











Fic. 2. Motion in a hyperbolic orbit between the Earth 
(A) and an inner planet (P). 
From which: 
n+q 
q 
In the same manner we can derive relations for outer planets. 
For ellipses: 
_ n—¢ 
colt 
q q 
For hyperbolae: 
q+1 
q 





(11) 


(12) 


(13) 


We can now set up a diagram (Fig. 3) with abscisa g and ordinate « for a 
certain n, say a voyage from the Earth to Venus. It is of advantage to use 
the logarithm of « because by so doing the more important area of the elliptic 
orbits is expanded and the less important area of the hyperbolic orbits is 


compressed. 


In this diagram every point inside the lines GH, HJ, KL and MN represents 


a possible orbit for a given . 
Fig. 3, 4, 5 and 6 are the diagrams for these voyages: 
from Earth to Venus: Fig. 3 
from Marsto Mercury: Fig. 4 
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hen San toMars: Fig.5 
from Earth to Mars: Fig. 6 


In the equations for the lines KL and HJ for outer planets the factor n does 
not appear. Therefore these lines are similar for all outer planets, although the 
length of the line HJ is different for different planets. The same rule applies 
to the line GH for inner planets. The values of m for different voyages are 
given in Table I. 


TABLE I,—Values of n for Various Voyages 

















n for voyages from 

Destination | Mercury Venus | Earth | Mars Jupiter | Saturn 
Mercury... | — 0-535 | 0-387 | 0-254 0-0744 | 0-0406 
Venus : 1-868 | — | 0-723 | 0-475 0-1390 | 0-0758 
Earth |} 2583 | 1-382 — 0-656 | 0-1923 | 0-1048 
Mars .. | 3-94 | 2-106 1-524 — 0-293 | 01597 
Jupiter | 1344 | 719 | 520 | 3-42 — | 0-545 
Saturn .. | 24-64 | 13:19 | os | O88 i] = 














Now, if we could express the values of the required energy, the duration 
of the voyage and the configuration of the planets at the moment of departure 
in terms of g, €, and ” 
only, we could calculate 
these values for every 
point in the diagrams 
and we would thus get 
a clear idea of the rather 
complicated circum- 
stances of interplanetary 
flight. 


4. The Orbit-to- 
Orbit Energy of 
a Voyage 








In Fig. 7 the velocities 
involved in the transfer 
of a spaceship from the 
orbit of planet A to the 
orbit of planet P are 


tial Fic. 7. Velocities involved in transfer from orbit of 
represented. planet A to orbit of planet P. 


These velocities are: 
V, the linear velocity of A in its circular orbit 
V, the linear velocity of P in its circular orbit. 
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V, the linear velocity of the spaceship in its orbit at point A. 
V, the linear velocity of the spaceship in its orbit at point P 


We know that: 


Vea .. ae es > + - 14 
r, (14) 
Vi2=° ., nee <3 on es ai 15 
p rs ( ) 


- 2 1 ’ 
vi=n(-—;) Ce aaa 
1 





Introducing: 
n="? and g= % 
r, nr 
we get: 
73_ Ve 
f= — Ks i wi on .. (18) 
%% — 
V2=V2. 41— i " im .. (19) 
% — 
V2=V.2. — o wae |, ee ae 


In order to adapt the initial velocity V, of the spaceship at the moment of 
departure to the required velocity V, in its interplanetary orbit we have to 
apply an extra velocity V ,. 

The direction of the velocity V, is perpendicular to the radius vector SA, 
the direction of the velocity V, is tangential to the orbit of the spaceship at A. 

The angle between the radius vector and the tangent is 0. 


Now we know that: 


dd 
tan@=—r— ae Ke “a ra 21 
an ra (21) 
And as: 
a(l — é?) 
= os ¥ ns - 22 
1+ e€cos¢ (22) 
We get: 


x 








r af al—e) ]_ ae(l — e)sing 
47 aliesces|- (1 + € cos ¢)? 
_ resing 
~ T+ ecos¢ uF a bee (23) 
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Thus 
yo .— + € Cos > =—tan@ .. ae ide a (24) 
dr € sin? 
Now: 
Via = V.2 + V;? — 2V,.V;. cos (0, — 32) = a (25) 
cos (0, — 4m) = sin 0, 
l > 
- ——__——— 26 
V (cot?@, + 1) ) 
Thus: 
2 =-{ 
cos (0, — 4m) = _ oa, +] 
“ (1 + € cos ¢,)? 
1 + 
+ <csg, _ a. 


VL + €? + 2e cos ¢,) 


Applying eq. (22) we find: 


1 + €cos v ated 
6ueniniain = , = Yi 2a(1 = €?) ‘ 
V (1 + & + 2 cos ¢,) FE ee ae Kee 
" 
_. B= 


Substituting eq. (28) and eq. (19) into eq. (25), we find: 


9 = a 
V ia re _ Ve . “4 a5 _ od i . 
q 


Vue ae See Lent —*) 
V q q 
; 1 , 
ees as “ — 
q 


l 
3—2 gl—é —- " ro 31 
alt M sf ” 


In order to aks the velocity V, of the spaceship at point P, to the velocity 
V, of the planet P we have to apply an extra velocity Vq,. 


l—ée 
— -- 29 
4 (29) 














or: 
Thus: 


In the same manner as above we can derive that: 


Vey _ fea ev i= sr _*} ey 
n 
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A measure for the energy required for the orbit-to-orbit transfer of a space- 
ship is the amount of fuel needed per unit of mass of the spaceship. 


This is expressed by the mass ratio R: 
M 
R= — - - - - ee 33 
M—m, ee) 
where M is the total mass of the ship and m, is the mass of the propellants. 
Now R is determined by the characteristic velocity V,,,. and the velocity 
of the exhaust-gases c: 





Ree So ook 


In this manner we can express the required energy as a function of the 
characteristic velocity, which is the sum of the separate velocities which the 
motor must confer to the spaceship. 


In our case: 


V en. = Vie i Ve» (35) 
Now we call: 
oe =E ss ”- - i .. (36) 


Thus E is a dimensionless number, which is a measure for the energy needed 
for the transfer of a spaceship from the orbit of a planet to the orbit of another 
planet. 

Then we find that for an elliptical orbit: 


/ — — €)/n) 
E=,/{s—2vigt—)- ‘} + {fen ~i} 
q n q 


Likewise we find for a hyperbolical orbit: 


E=,/{3—2viae ! +i bs, /{ t= eve Mit) @g 
q n q } 


These equations are valid for inner as well as for outer planets. 

In this manner we have expressed the required energy of an orbit in terms of 
q, « and only. 

We can now calculate the values of E for several points in the above- 
mentioned diagrams and then trace lines of equal energy in these diagrams. 

We call these lines iso-ergs. Thus an iso-erg gives a survey of all possible 
orbits with the same energy-requirement. 

Fig. 8 shows the iso-ergs for voyages from the Earth to Venus. We see 
that E rapidly increases for ellipses of increasing eccentricity and more slowly 
for increasing g. For hyperbolic orbits the circumstances are more complicated. 

Furthermore we see that some elliptical orbits of high eccentricity require 
more energy than some hyperbolic orbits of low eccentricity.* 

* Though the author calculated about one hundred points in the diagrams, this number 


is still too small to trace the lines accurately. Therefore the lines must be considered as 
showing only the general trend of the parameters. 
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ISO-ERG DIAGRAM 
EARTH—> VENUS 
ROUTE A 





1-0 12 1-4 
q 


Fic. 8. Iso-erg diagram for voyages from Earth to Venus (Route A). 


When interplanetary flight is fully developed there will doubtless be a place 
for slow as well as for fast transport, even as in the present time slow transport 
by ship exists beside fast transport by airplane. 

We can imagine big freighters, which start at a date when the configuration 
of the planets is favourable for a low-energy-orbit, and which drift slowly 
along their low-eccentric orbit to the planet of their destination. On the other 
hand small ships on urgent business must be able to start at any moment and 
reach their destination in a short time at the cost of a multiple of the energy 
of the former ships. 

The total energy for an interplanetary voyage is much higher than the 
energy needed for the transfer from orbit to orbit. We must add: 

(1) The energy needed to escape from the planet of departure. Although 
the escape velocity is in reality the velocity necessary to escape to 
infinity the difference for interplanetary distances is small. 

(2) The energy needed to brake the fall on the planet of arrival. This 
energy is also practically equal to the escape velocity of the planet. 

(3) The extra velocity needed to compensate for the fact that the spaceship 
starts from the planet with zero velocity with a finite acceleration. 

(4) A similar extra velocity must be added to the escape velocity of the 
planet of arrival. 
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(5) In the case of a start from a planet with an atmosphere the loss of 
velocity, owing to friction of the atmosphere on the spaceship must be 
taken into account. 
Taking as an example a voyage from the surface of the Earth to Venus we 
get the following figures for an orbit with E = 1-0: 


Velocity for the transfer from orbit to orbit 


(1-0 x 29-8 km. /sec.) so af as 29- S km. /sec. 
. Escape velocity of the Earth + in 11-2 km. /sec 
2. Escape velocity of Venus... 10-4 km. /sec. 
3. Extra velocity for a mean acceleration of 5g at 
start oe 1-3 km. /sec 
4. Extra velocity ‘for a mean ‘deceleration of 5g at 
arrival me ee ca 4 +2 km. /sec. 
5. Air resistance in the atmosphere of the Earth 5 km. /sec. 


Total characteristic velocity a _ ws 54-4 km. /sec 


We see that the energy for the transfer from orbit to orbit is only 55 per 
cent. of the total energy. 

In this calculation we did not take into account the possibility of braking 
our fall in the atmosphere of Venus, because this possibility is still a moot point. 

If the voyage is from a spacestation circling the Earth to a spacestation 
circling Venus, the required energy is of course much less, but in this case we 
have to add the energy for the transfer from the surface of the planet to the 
spacestation and vice versa. 


5. The Duration of the Voyage 
The general equation for a curvilinear movement in polar co-ordinates is: 


v2=rr(g?+ (2 .. ol a .. (39) 


For an elliptical orbit this equation becomes: 


= -u( > — ;) éa ‘i ™ .. (40) 
of ') 7 yr .. (41) 
L a 


From which: 





rr a 
Furthermore: 
sd a(l — €*) 
"Fie cos d 
From which: 
a een 
é= emo" i ae .. (42) 
ve 
and 
b= — oVG. — 3.9 J .. (43) 


r4/(2ar — r? — a(1 — )) 
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Equation (43) substituted in eq. (41) gives after simplification : 
a= (*). aaa - Sug 
\u 


V/{ — r? + 2ar — a (1 — &)] 
When the object is at perihelion the time ¢ = 0. Its radius vector is then: 
yo = a(l —e). 


After integration we then find that the time ¢ for the object to arrive from 
a point of the orbit with radius vector 7 to the perihelion is: 


t= a’ > pt { arc Cos (*=*) — ll _ (‘ ~ “y'|} -» (45) 


We call: 
arc COs (3 ") =F .. mh ee ve (46) 
ae 


t a®. u-i(r — esinz) as ia ne (47) 





Then: 


The duration for a voyage from a point on the orbit with radius vector 7,, to 
a point with radius vector 7, is: 
3 
T=t,-—t,=a?. p[(7, — 72) — e(sint, —sin7,)] .. .. (48) 


where: 
a—r —1l 
T, = arc cos a) = are cos ( / ) me .. (49) 
ae ge 


T; = arc cos (F — *) = arc cos (*—") m -. (50) 
ae ge 


The time for a complete revolution of the Earth around the Sun is: 





and: 








T, =79,'. p+. 2a ws i i“ a” = 
Dividing eq. (48) by eq. (51) we find the duration of the voyage in years: 
T = a [(r; — tT.) — e(Sint, — sin7,)] (years) .. és -- (52) 
or: 
T = 58-13 q [((r, — tT.) — €(sin rt, — sin 7,)] (days) is -- (53) 
In the same manner we can derive for a hyperbolic orbit: 
T = 5813 g3[¢(sinh 7, — sinh 7,) — (r, — 7)] (days) (64) 
where: 
™1= are cosh 4 ; eh va -. (55) 
a T. = arc cosh i= hy - .. (56) 


ge 
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Fic. 9. Isochrone diagram for voyages from Earth to Venus 
(Route A). 


In this manner we have expressed the duration of the voyage in terms of gq, 
«and only. In our diagram we can now draw lines of equal duration, which 
we Call isochrones. 

In Fig. 9 are shown the isochrones for voyages from the Earth to Venus in 
days. The longest duration is 146 days, but by using high-energy orbits this 
time can be reduced to some days. 

If we use a diagram which contains isochrones as well as iso-ergs we can 
easily solve problems such as the following: 


“Find the orbit which requires the least energy for a given duration.” 





Examples. 
= (days) Ewin q , 
40 0-78 0-99 0-38 
50 0-61 0-94 0-32 ) 
60 0-49 0-92 0-27 


Or, conversely: 
“Find the fastest orbit for a given energy.” 


Examples. 


E T min. (days) q € 

0-8 39 1-00 0-39 
0-7 44 0-97 0-36 
0-6 51 0-94 0-32 


0-5 59 
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These problems would entail a considerable quantity of computing work 
to solve them, the more as equation (52) is transcendent, so that they cannot 
be solved algebraically. 


6. The Configuration of the Planets at the Time of Departure 

As the planet of our destination moves around the sun with an angular 
velocity which is different from the angular velocity of the planet of departure, 
the configuration of the two planets changes continually. 

The configuration at the time of departure determines in part the orbit 
which we have to follow in order to reach the planet at the exact time and 
place. 

If the radius of the orbit of the planet P is m in A.U. and the radius of the 
orbit of the planet A is m in A.U., then P will make one revolution around the 


Sun in n2 years and A in m: years. 
—3 —3 
In one year P describes an angle of 27m 2 radians and A an angle of 27.m 2 
radians. 
Then the variation of the angle between the two planets will be per year: 


+ 22(n-*/* — m-*) radians... a a> 


When the angle between the two planets has increased with 27 radians 
we have the same configuration as at the beginning. 
The time in years between two similar configurations is thus: 


. 1 
T°? = + —»_——_, years... “s -» (58) 
n—*!? — m-*?2 


When 4 is the Earth, the time between two similar configurations is: 
l 


es 
n~*i= — J 


7? = years .. - +“ .. (59) 


(+ for inner planets, — for outer planets.) 


For the Earth and one of the following planets the times 7* are given in Table II. 








TABLE II 
Planet n T* (days) | Variation of ¥ per day 
Mercury " ool 0-387 116 — 3-10° 
Venus as os] 0-723 | 584 — 0-617° 
Mars .. - oat 1-524 780 + 0-462° 
Jupiter a « 5-20 398 + 0-904° 
Saturn . on 9-54 379 + 0-949° 














The times T* are of course the so-called synodic periods of the planets. 

We see that the waiting times for a similar configuration are longest for the 
nearer planets. 

The calculation of the configuration of the planets goes as follows: 

In Fig. 1 the configuration of the planets is represented by the angle ¥. 
This angle can be positive or negative. We call ¥ positive when P, is situated 














342 M. VERTREGT 


before A if one goes in an anti-clockwise direction from the perihelion of the 
orbit of the spaceship (the planets move around the Sun in an anti-clockwise 
direction). 
Furthermore: 
Y- wa) Pi 180°): 

When T is the time needed for the spaceship to travel from A to P,, and T, 
is the time needed for the planet P to go from P, to P,, then obviously T must 
be equal to T, for the spaceship to reach the planet at the exact time and place. 


Now: 
3/2 
T=T,= > (@+ P) (years) _ .. (60) 
a7 
From which: 
Y = 2nn-*” . T — @ (radians) .. “se .. (61) 
where: 
_ - — 
® = $, — $2 = arc cos [e—] — arc cos (“2 "| (62) 
€ ne 


As T is known from eq. (52) we have expressed Y in g, e and” only. We 
can now trace lines of equal configuration in our diagram, which we call tsogons 
(Fig. 10). For g < 2-0 and for elliptical orbits the isogons range from about 
+ 50° to about 0° for a voyage from Earth to Venus. 
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Fic. 10. Isogon diagram for voyages from Earth to Venus 
(Route A). 
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Now, if we trace the iso-ergs, isochrones and isogons in one diagram on a 
sufficiently large scale, using different colours, say red, green and black for 
greater clearness, we can easily solve many intricate problems of interplanetary 
travel, which will certainly arise when spaceflight is a reality. 

For instance: 

Example 1. 

A spaceship wishes to leave Earth for Venus at a certain date when ¥ =+ 14°. It 
is the intention to arrive at Venus in 40 days. Which E is needed for this voyage and 
which are the elements of the orbit? 

By following the isogon of + 14° we meet the isochrone of 40 days at a point with the 
co-ordinates: 

q = 0-984 and ¢« = 0-385. 
The E for these co-ordinates is 0-78. 


Example 2. 
At 15 May, 2058 the ’ of the Earth and Venus is + 90°. A spaceship with a maximum 
E of 1-0 wishes to depart at the earliest opportunity. When can the spaceship start and 
when will it reach Venus? 
As ¥ diminishes 0-617° per day (see Table II) the earliest date is determined by the 
highest possible ’. 
By following the iso-erg E = 1-0 we find that the highest ¥’, situated at the extreme 
left of the diagram is + 45°. 
Thus it will take: 
90 — 45 
——_.- days = 73 days 
— 
before the spaceship can leave. The co-ordinates are: 
q = 0-640 and « = 0-562. 
Looking at the isochrones we find that 
T = 57 days. 


Thus the ship can arrive at Venus 73 + 57 = 130 days after 15 May, or at 22 September. 


7. Hohmann Orbits 

These orbits, called after the German pioneer of spaceflight, Dr.-Ing. 
Walter Hohmann who first proposed them, deserve special attention, because 
they will probably be used in the first endeavours to reach the planets. 

The Hohmann-orbit is represented in our diagrams by the lowest point H 
(Fig. 3). We see in the iso-erg diagram (Fig. 8), and it can be proved mathe- 
matically, that this orbit requires the minimum energy, but at the same time 
the isochrone-diagram (Fig. 9) shows that the duration of the voyage is longest. 

The Hohmann-orbits are co-tangential to the orbit of the planet of departure 
as well as to the orbit of the planet of arrival. Therefore, assuming that the 
orbits of the planets are circular the angle ® of these orbits is 180° and for 
inner planets ¢, = 180° and ¢, = 0, whereas for outer planets ¢, = 0° and 
¢. = 180°. 

The elements of the Hohmann-orbits can be easily calculated from eq. (6). 

We have in this case for inner planets: 
sor Fo Te ae 


€e= 





q 
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from which: 








a4 as 
¢= = and ply. an i se .. (64) 
2 lin 
For outer planets: 
n+ 1 n—l1 - 
q - and €= — _ oe .. (65) 
2 n+ 1 


Substituting there values in eq. (36) we get for inner <ayaee 


“ “ 
Eun ees ll I+ Rie: a )-1 .. (66) 
n -+ Vn nm + 7 


For outer planets: 


E = VG 7" n I+ val VG) (67) 
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Fic. 11. Measure of the energy required for Hohmann orbits to various planets. 


Fig. 11 shows E for different values of m. E has a maximum for » ~ 15 
and then slowly decreases to 1/2 — 1 = 0-414 for n > oo. 

An object describing an ellipse with semi-major axis g around the Sun 
makes a complete revolution in g*/? years. 

As a spaceship following a Hohmann-orbit makes exactly half a revolution, 
the duration of the voyage follows from eq. (64) and (65): 


T=3("5 an - 8 (years) - ei a 


For the planet P we have according to eq. (60): 
3/2 
-— (® + W) (years) ¥ Te 





As: 


@=n 
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we get from eq. (68) and (60): 
Y= “| : -) ie | (radians) .. — 





2n 
Table III gives some particulars about the Hohmann-orbits to the planets: 


TaBLeE III.—Data for Hohmann-type Orbits 











| 1 

From Earth to: n q € E years days Y 

Mercury -- | 0387 0-694 | O441 | 0-575 | 0-289 105 | + 252 
Venus .. .. | 0723 | 0-862 | 0-161 | 0-175 | 0-400 146 + 543 
Mars .. ¥ 1-524 1-262 | 0-208 0-187 | 0-709 259 — 443 
Jupiter _ 5-20 37102 | 0-678 | 0-485 | 2-732 — — 971 
Saturn .. | 9°54 527 | 0810 | 0-528 | 6-05 _ — 106-0 
Uranus -- | 19°19 10-10 0-901 | 0-535 | 16-12 - -1l1-4 
Neptune .. | 30-07 15-54 | 0-936 0-527 | 30-6 — 113-2 
Pluto .. ++ | 395 20-2 | 0-950 | 0-520 | 45-5 - —114-0 





These figures show that the duration of a voyage to a far planet is excessively 
long. If Pluto will ever be visited by human beings they will have arrived 
there along a faster orbit. 

Though the Hohmann-orbits are the most economic from a point of fuel 
consumption they are not the easiest to navigate. It is clear that in the case of 
a short high-energy orbit a small deviation at the moment of departure from 
the calculated magnitude and the direction of the thrust will give only a small 
error in place and time of arrival. 


But equally small deviations in the case of the long and slow Hohmann- 
orbits can cause that the spaceship misses the planet of destination altogether. 


Furthermore we must remember that E represents only a part of the total 
energy required for the voyage. As in the case of the other orbits we have 
to add the escape-velocities, etc., to the orbit-to-orbit velocity. 

In Section 4 we calculated these extra velocities for a voyage from Earth to 
Venus at 24-6 km. /sec. 

Table IV shows the ratio of the total energy of a high-energy orbit to the 
total energy of a Hohmann-orbit for the case Earth-Venus for different values 


of E. 








TABLE IV 
E Ratio 
0-175 1-0 (Hohmann-orbit) 
0-5 | 1-32 
1-0 | 1-83 
1-5 | 2-33 
2-0 | 2-83 
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Fic. 12. Alternative elliptical routes between the orbits of 
two planets. 


We see that the total energy for a high-energy orbit does not increase so 
much as would be expected from the orbit-to-orbit energy alone. For instance, 
an orbit with an E of ten times the E for the Hohmann-orbit requires a total 
energy of only two and a half times the total energy of the Hohmann-orbit. 


8. Indirect Routes 
Till now we have considered only the possibility of reaching a planet by the 

shortest route AP (Fig. 12), which we call Route A. We saw that the configura- 
tion ¥ for this route lies between + 50° and about 0°, if we follow orbits of 
reasonably low energies. This would mean that in the case of the Earth-Venus 
journey the latter planet would be accessible only during a period of about 
80 days out of 584 days. But as an elliptic orbit generally intersects the 
circular orbits of the two planets in four points we have three alternative 
routes, namely: 

Route B: A, B, P’ (Perihelion Route) 

Route C: A’, C, P (Aphelion Route) 

Route D: A’,C,B, P’ (Indirect Route) 


As point P’ is symmetrical to point P in relation to the major axis and as A’ 
is symmetrical to A, it is clear that the triangle of the velocities (Fig. 7) for 
point P’ will be exactly equal to the triangle of the velocities for point P, and 
those for point A’ equal to those for point A. Therefore the energies for these 
routes will be equal to the energy for Route A. 

The duration of the voyage will of course be longer than for Route A, the 
difference depending on angle ¢, for Route B and on ¢, for Route C. 


In connection with the difference of the duration of the voyages, the angle 
Y will also be different for each route. 
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When 7; is the time for the spaceship on Route B to cover the angular 
distance: 


®,= ©, + 24.=¢4,+% ‘. aa boa (71) 


and 7, is the time for the planet to cover the distance ®, + Wp, then: 


n3!2 n3/2 re 
Tz = Se = on (DO, a P,) — on (¢; — - db. — ¥,) ee ee (72) 

From which: 
Y, = 2r.n-*".T, — (d,+¢4) .«- a "ys (73) 


When 7, is the time for the spaceship on Route C to cover the distance: 


®. = ®, + 2(7 — ¢;) 


= 27 — (¢, +d.) .. % ‘4 oo ae 
and T, is the time for the planet to cover the distance ®, + Wz, then: 
3/2 3/2 
Tc = T, = — (©, + P.) = — [2 — (fr +a) + Pc) .. (78) 
«7 a7 
From which: 
Y, = 2n.n-*".T, — [2a — (6, + 4,)]_ -- .. (76) 
When 7p is the time for the spaceship to cover the distance: 
®, = 2a — D, = 2m — (d, — de) i _.. cor 
and T, is the time for the planet to cover the distance: 
®, + VY, = 27 — (¢, — 4.) + ¥, i xs. oe 
then: 
3/2 . n3/2 - | me 
T,=T,= Sa (®, + ¥,) = ._ we ($1 — $2) + Po) e (79) 
From which: 
Y, = 2n.n-*?.T, — [24 — (6, —4s)]_-- .. (80) 
It is easy to derive from eq. (52) and eq. (71), (74) and (77), that: 
3/2 
t,o 7,4 to, ete 260.) aoe ee 
7 
3/2 
T,=1,+ _ (7 — 7, + €sin7,) ke -- (82) 
T,=-—T,+¢" .. ok - in .. (83) 


(All angles in the foregoing equations are in radians, all times in years.) 


For hyperbolic orbits we have only one alternative route (Fig. 13), namely: 
Route B: ABP’ 
and we find that: 
/2 


T,=7,+£ (e sinh tT, — 7.) - .. (84) 


Y, = ten, ~-b+0b oe 
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Fic. 13. Alternative hyperbolic routes between the 
orbits of two planets. 
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Fic. 14. Isochrone diagram for voyages from Earth to Venus 


(Route B). 
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Fic. 15. Isogon diagram for voyages from Earth to Venus 
(Route B). 


We can now make diagrams for the isochrones and isogons of Route B 
(Fig. 14 and Fig. 15) and of Route C (Fig. 16 and Fig. 17). 

Route D is of little importance because of its long duration; the practical 
routes would be those of which @ is about 180°, thus the orbits in the vicinity 
of the Hohmann-point H (Fig. 3). 
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Fic. 16. Isochrone diagram for voyages from Earth to Venus 
(Route C). 
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Fic. 17. Isogon diagram for voyages from Earth to Venus 
(Route C). 


In the diagram of Route B a part is left blank, where the spaceship at 
perihelion would come too close to the Sun. We took for the minimum dis- 
tance about half the radius of the orbit of Mercury, viz.: n = 0-2. 

In the diagram of Route C that part is left blank where the duration of the 
voyage would exceed one year. 

It will be seen that the configurations Y for the Routes B and C comprise 
all angles from + 180° to — 180°. From this it follows that it is possible to 
start for a planet at any desired moment. 

This is important, because we saw earlier that using Route A the possible 
starting time is only during less than three months in a period of more than 
one and a half years. 

The advantage of an alternative route is demonstrated by the fact that in 
Example 2 (Section 6) we needed 130 days, including a waiting time of 73 days, 
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Fic. 18. Relative configuration of Earth and Venus for the years 2057 and 2058. 
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to arrive at Venus. By following Route C we can start immediately and 
arrive at Venus, using the same energy, after 100 days, that is 30 days earlier 
than by following Route A. 

Fig. 18 shows the angle Y for the case Earth-Venus during the years 
2057 and 2058. 

From the isogon-diagrams we can see that the favourable time for Route A 
is during the months January and February 2057, the favourable time for 
Route B is from March to September 2057 and the time for Route C is from 
October 2057 to July 2058, when Route A can again be followed. 


Of course, the periods for the different routes overlap in part and are 
furthermore to some extent dependant on E. 


As a result of the fact that the duration of the voyage for Route C is longer 
in the beginning of the period and becomes steadily shorter at a later date, the 
spaceships which start from Earth during this period will all arrive at Venus 
shortly after one another. 

The most unfavourable time to start from Earth is the beginning of the 
period for Route C, when the voyage lasts about eight months. 

The low-energy ships from Earth which follow Hohmann-orbits will arrive 
at Venus in May 2057 and December 2058, in general when ¥ =— 36°, after 
a voyage of 146 days. 


9. Conclusion 

In the introduction to this article we made several assumptions in order 
to simplify the calculations. One of these assumptions was that the planets 
describe circles around the Sun. The eccentricities of the orbits of the planets 
are as shown in Table V. 








TABLE V 
Eccentricity | Inclination of the orbit 
Planet | of the orbit to the plane of the ecliptic 
| 
Mercury .. ja ¥ 0-206 7-0 
Venus... ae 0-007 3-4 
Earth 0-017 0 
Mars 0-093 1-9° 
Jupiter — 0-048 1-3 
Saturn... ie J 0-056 2-5 
Uranus 0-047 0-8 
Neptune .. 0-009 1-8° 
Pluto 0-248 17-2° 











Except for Mercury and Pluto and to a lesser degree for Mars, the deviations 
from a true circle are small, and the resulting errors for the case of the Earth 
and Venus are negligible in approximate calculations. 


This is also true for the second assumption that the planes of the orbits of 
the planets coincide with the plane of the ecliptic. 
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As Table V shows, the inclinations are small except in the case of Mercury 
and Pluto. 

It is possible to design the above-mentioned diagrams based on the real 
circumstances, including the perturbations caused by the attraction of the 
planets, but I do not think that the enormous amount of computing work 
necessary to introduce these corrections would be worth while, the more as 
this work should be done anew every year. The advantage of the approximate 
diagrams is that they can be calculated once and for all. 

The calculation of the exact path of a spaceship is so complicated that this 
work must be done for every voyage before departure with electronic com- 
puting machines. In the future special computing machines will doubtlessly be 
constructed for this purpose and the accurate computing of the elements of an 
orbit: the exact time of departure, the exact magnitude and direction of the 
thrust will be routine work for these machines and will occupy only a short time. 

The navigator of the spaceship must have at his disposal the means to 
compare the real orbit with the calculated orbit and to apply small corrections 
if necessary. Special techniques must be worked out so that these corrections 
can be computed on board the ship with simple calculating apparatus. 

The author does not suppose that the diagrams, developed in this article, 
can supplant these intricate calculations. No diagram could be designed 
exactly enough to satisfy the required accuracy. But the practical utility of 
the diagrams is provided by the clear survey they give of the possibilities and 
particulars of interplanetary voyages. They will facilitate the provisional 
planning of a voyage for given circumstances, they show the possible alterna- 
tives and they give approximate answers to questions as to the favourable 
configurations of the planets, the required energy and the duration of a voyage. 
For this special purpose the diagrams will prove their usefulness in inter- 
planetary flight. 

REFERENCE 
(1) M. Vertregt, J.B.I.S., 1956, 15, 234. 


Addendum+ 
Some time after concluding this paper the author discovered that a con- 
siderable simplification can be affected by not using g = “ for the abscissae in 
"1 
the diagrams, but by using: 
a(1—e?) 


" 





a(1l — e*) is the “parameter” of a conic section; it is the radius vector with a 
perihelion-angle of 90°. 
The advantages of using / are: 
(1) The boundary-lines GH and HJ—/L of the diagrams become straight 
lines, if the diagram is designed on ordinary squared paper (see Fig. 19 


t Received 30 October, 1957. 





(2) 
(3) 


(4) 


E=-,/[s—2ve-(5# 
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for Earth-Venus and Earth-Mars voyages). The angles of these lines 


with the ordinate are: 


angle inner planets outer planets 
7 45° arc tan 
B arc tan 45° 


All possible elliptical orbits are comprised in the triangle GH J. 


The parabolae which did not fit into the former diagrams, now fit in 
nicely. They are represented by the line GJ. The length of this line 
is 2m for inner planets and 2 for outer planets. 


The equation for the orbit-to-orbit energy now becomes: 


+ {fA weet 
bal” Hasan ees ac 


and is valid for all conic sections. For parabolae: « = 1-0. 
As q has disappeared from this equation, there is no discontinuity 
for g = 00 (parabolae). 


As a result of this the curves of the diagrams now run without interrup- 
tion from the area of the ellipses, through the line of the parabolae to 
the area of the hyperbolae. 
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As an example the new iso-erg diagram for Earth-Venus is shown 
in Fig. 20. A comparison with Fig. 8 shows the advantages of this 


method. 


ERRATA 


Orientation in Space 


The following corrections should be made to the paper by M. Vertregt 
(J.B.I.S., November—December, 1956, 15 (6), 324-338). 

: = ——, send ¢ = : — 
—cos B’, cos B’;’ ~ cos(— B’;) cos B’; 


Page 332, left hand side of equation (15). For tan 8, read tan 6. 





Page 327, line 8. For-:r - 


Page 333, line 13 from bottom. Fortan=6@...,readtanOd=.... 
Page 334, line 7. For SA + Q — 95, read SALQ—DO. 

Page 334, line 8. For AR + SA,read ARILSA. 

Page 336, line 3. For @ = 0, readd = 0. 
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A SIMPLE METHOD OF PLOTTING THE 
TRACK OF AN EARTH SATELLITE* 


By R. A. ANDERSON,t PH.D. AND C. S. L. Keay,t M.Sc. 


SUMMARY 


A technique is described whereby interested persons may plot from limited information 
the track of a satellite in an orbit of low ellipticity and make short term predictions without 
recourse to complicated computing procedures. Simplicity and speed of application have 
been aimed at, rather than a high degree of accuracy. 


1. Introduction 

When the announcement was made of the launching of the first Russian 
satellite no reliable information about its path or times of appearance over 
New Zealand was available. A simple method of plotting the track of the 
satellite was therefore devised, enabling predictions to be made of the times 
at which a radio and visual watch should be kept. These predictions resulted 
in many more visual observations than would otherwise have been made. 
The importance of these observations was emphasised when it became apparent 
that the visual object was, in fact, the third-stage rocket, which had no radio. 


It is the purpose of this paper to describe the method of plotting the tracks 
of satellite objects in order that interested persons may make their own pre- 
dictions. A further discussion will indicate how local sightings resulting from 
these predictions may be used to derive more accurate track positions, thus 
improving the accuracy of later predictions. 


Such a project would be well within the capabilities of amateur groups such 
as small astronomical societies or advanced secondary school classes. 


2. The Basis of the Method 

The method consists of preparing a template which represents the approxi- 
mate path of the orbit (assumed to be circular) projected on the Earth’s surface. 
By means of this, any reports (or authoritative predictions) of the satellite 
position may be related to the time and longitude at which it crossed the 
Equator (the node of the orbit). From these times an ephemeris may be 
compiled. 

The information required to construct the template is the period, 7, and 
inclination, 7, of the orbit, which, if unknown, may be derived from at least 
two (but preferably more) reliable reports of the satellite’s position. 


The simplest method of deriving these orbital parameters is to obtain first 


* Manuscript received 25 November, 1957. 

+ Formerly at University of Canterbury, Christchurch, New Zealand; now at University 
of Western Australia, Nedlands, Western Australia. 

¢t University of Canterbury, Christchurch, New Zealand. 
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the angular distance around the orbit, A@ between the two reported positions. 
For this we use the relation 





cos A@ = sin ¢, sin gd, + cos¢, cosg,cos{ .. ie “ (1) 
where 
¢,,¢, = latitudes of first and second reported positions 
t—t 
C = ),—A,+- z * degrees 
A,, A, = longitudes of first and second reported positions. 


t,, tg = first and second reported times in minutes (U.T.). The last term 
in the equation for { represents the correction due to the rotation of the Earth 
and the sign is taken positive for direct orbits. 

In equation (1), A@ must always be taken in the same quadrant as ¢. 

Once A@ has been derived the inclination and period follow immediately 
from equations (2) and (3): 
cos , Cos dy sin [ 





1 9 
COS 2 sin AQ (2) 
(if A@ lies within 10° of 0° or 180°, an unreliable value of ¢ may result) 
tg —t 
T =. ——__... ~ we «| = 


In equation (3), » represents the integral number of revolutions of the 
satellite between the reported positions. Generally this will be unknown and 
has to be chosen to give sensible and consistent values of T. 

To illustrate the procedure typical results obtained for the first Soviet 
Sputnik are given in Table I. However, the values of i and T were available 
for Sputnik I at an early stage, which considerably simplified the problem. 


TABLE I.—Results Derived from Predicted Positions of the First Russian Satellite 





























Times Longitude | T 
Date | Places (U.T.) and latitude r4 Aé i |min.| n 
Jp | — satel Babe Coca UES! sesdets 
6 Oct. | 1. Baghdad 19 hr. 34min. | 44° E, 33° N ° o| aro | 9 
1957 | 2. Hiroshima | 23hr. 13min. | 132° E, 35° N | '49°| 104°] 65° | 96 | 2 
Sey gee erty pooragp eer: aa: ae | an, wh ae 
8 Oct. | 1. New York | Ilhr. 59min. | 287° E,41°N 20 | moo | gro | gg | 
1957 | 2. Riode | 12hr. 18min. | 318°E,23°s | 36° | 72° | 65° | 94 | 0 
Janiero | 
09 hr. 38 min. | 145° E, 39°S 216° | 260° | 67° 9% | 8 




















§ aS Cee 
8 Oct. | 1. Melbourne | 
| 2. Sydney | 23hr. 35min. | 152° E, 34°S 





With this information a template may be constructed as follows. First 
a table is derived giving values of ¢, the satellite’s true latitude referred to the 
Earth’s equator, and 2’, the satellite’s longitude referred to the node of the 
orbit projected on a non-rotating Earth in terms of 0, the angular distance 
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along the orbit from the node. Since the orbit is assumed circular we may 
write 6 = 360t/T. Fig. 1 shows the relation between 6, ¢ and 2’. 


This may be expressed in the following equations 
sing@=sinisin@ .. a ‘ “ (4) 
tan A’ = cos i tan 6 _ vs ‘a (E 


i) | 
~— 


It is now necessary to include a correction for the rotation of the Earth, 
yielding true values for the satellite’s instantaneous longitude: 
6T 


aa. J a — 


A=N 
where TJ is expressed in minutes and @ in degrees. The negative sign must be 
taken for direct orbits. The factor 26 is included to account for nodal regression 
but is neglected in deriving the template. When the values of ¢ and A have 
been tabulated as a function of ¢ up to 7/2, the outline of the template may be 
plotted, and a time scale 
marked around the edge. 
The size and shape of the 
template depends on the 
map projection, which 
must have symmetry 
with respect to the 
Earth’s axis, in order 
that the template will be 
accurate at any longitude. 
A Mercator’s projection 
is suitable and for this 
the template will approxi- 
mate to a half cycle of 





Fic. 1. The orbital path of a satellite as projected on a ’ =". 
non-rotating Earth. a sine wave. This is 
illustrated in Fig. 2. 


3. Application of the Method 

An ephemeris may now be derived by fitting the template to reported 
positions of the satellite and noting the times and longitudes of nodal crossings. 
In order to identify correctly whether the nodal crossing is ascending or descend- 
ing, it is better to commence if possible with a pair of positions which must 
obviously lie on the same revolution (e.g., New York and Rio de Janiero in 
Table I). All other positions may be related to nodal crossings; care being 
taken to ensure that the times at the nodes differ by mT for like nodes and 
(nm + 4)T for unlike nodes. When the difference in time and longitude between 
successive nodal crossings has been determined in this manner to a reasonable 
degree of accuracy it is possible to derive tracks and times of appearance over 
any part of the world and make predictions for a day or two in advance valid 
enough to enable watches, both radio and visual, to be kept. Periodic correc- 
tions, preferably made from local observations, will of course, be necessary 
as time goes on. 
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4. Utilisation of Local Observations 

Whenever a satellite proves to be visible to the naked eye, the most 
useful type of observation which can be made without special equipment 
is to note accurately the times of transit of the satellite vertically above 
or below convenient stars. This does require a reasonable knowledge of 
the main constellations and the use of a star map to ensure positive identifica- 
tion of the stars used. If at least three such observations can be made each 
time the satellite appears it is possible with the aid of a nautical almanac and 
standard astronomical procedures to plot its track. Unless these observations 
can be made to a high degree of accuracy it is better to use as far as possible 
graphical methods, which will now be described. 

An irreducible minimum of computation consists in finding the local 
azimuth of the check stars at the time of observation, as shown in Appendix I. 
These azimuth bearings are then drawn on a map covering an area extending 
up to about 1,000 miles around the observation post. The best type of map 
to use is an orthomorphic conic projection on which the azimuth bearings will 
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Fic. 3. Diagram showing how the projected tracks are fitted to times of transit 
T,, T, above Venus and 7;, 7, below a Aquila as observed from two widely spaced 
stations. In this case the time differences 7,—7,, T,—T; yield the track speed of 
the satellite and the time differences T,—7,, T7,—T7,, T,—T,, Ty—T,, when converted 
to track lengths, are used to locate the actual track position. 
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closely approximate to straight lines, rather than rhumb lines as with a Mercator 
projection (see Fig. 3). 

A new template and time scale, covering the observed portion of the orbit, 
is required for drawing in the tracks on the local map. Although the map 
projection may be different, the method described in Section 2 may still be 
used. For an orthomorphic projection the outline of the template will approxi- 
mate to a straight line, which may then be moved along the latitude lines 
until the track differences fit the azimuth bearings already plotted. The 
track differences may be easily obtained from the time differences between 
observations if the track speed of the satellite is known: this may be found to 
a first approximation by dividing the Earth’s circumference by the period of 
revolution of the satellite. If the tracks obtained for different pairs of times 
(and bearings) do not line up with one another a modified value of the orbital 
inclination may be indicated. If so, an improved template may be constructed. 

When at least two complete tracks have been derived it is possible to 
correct the ephemeris and thus improve the accuracy of future predictions. 
Too much must not be expected of this method, however, particularly if the 
orbit is fairly elliptical. It may be noted that this means of deriving an 
ephemeris automatically includes such an effect as a regression of the nodes. 
Similarly, changes in the period of the orbit will make little difference to the 
form of the template and track directions but must be considered in determining 
the nodal positions and times. 

Using the methods described above, the authors derived tracks in the 
vicinity of New Zealand from visual observations made on the evenings of 
October 9 and 12, 1957. These have been reported elsewhere,’ and it is interest- 
ing to note that the regression of the nodes was found to be 11’ 36” + 48” per 
revolution. The regression expected due to the oblateness of the Earth may 
be derived from the theory of Spitzer,? and of Blitzer, Weisfeld and Wheelon,*4 
and is found to be 12’ 41” + 30”. Thus it appears that there may be a signifi- 
cant difference between the observed and theoretical values. 
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Appendix I 


The steps in computing the local azimuth of a star from its right ascension 
and declination at a particular time are given here for convenience. The 
symbols used are: 


= right ascension. 

= declination. 

= longitude of observing station (East). 
= latitude of observing station. 

= azimuth. 

= equation of time. 

= hour angle. 

= zenith distance (90° — altitude). 
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1. Q&star — %sun = Aa. 

2. Local transit of sun = 12 hr. — A — E (U.T.) (where A and E are 
expressed in hours). 

3. Local transit of star = local transit of sun + Aa« (U.T.) (again, Ax must 
be in hours). 


4. Hour angle of star, h = 
local transit of star (U.T.) — local time of observation (U.T.). 


5. cosz=sin Ssind + cos 8 cos¢ cosh. 


f cos 6 sin h. 
6. sin A = ——— 
sin z 
For negative values of h, A is West of North. 
For positive values of 4, A is East of North. 
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THE SATELLITE TELESCOPE* 
By F. A. Smitu,+ Grad.I.Mech.E., FELLow 


SUMMARY 


Some of the factors affecting the design of satellite telescopes are discussed, and a 
suggested layout for such a telescope is put forward. 


Introduction 

When satellite stations have been set up in orbit around the Earth, one 
of the first tasks to be undertaken will be the construction of an extra-terrestrial 
observatory. The advantages of such an observatory have been set out by 
C. A. Cross! The question of its design has attracted the attention of a 
number of authors, but very often their proposals have their foundations in 
some ingenious but otherwise uncertain principle. Whilst such ideas may 
become practicable in time, it is felt that the first observatory will be made 
with as few untried features as possible. Such an observatory would be set 
apart from the satellite station in order to be free from any rotational 
movements or vibrations. 

The first major decision to be faced is whether the telescope is to be a 
reflector or a refractor. The more protracted grinding and figuring of the 
latter, and the chromatic perfection of the former tilt the scales heavily in 


* Manuscript received 23 May, 1957. 
+ Bristol Aero Engines, Ltd. 
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favour of the reflector. To provide a similar degree of meteor screening to 
a refractor as can be more readily supplied to a reflector gives rise to an 
appreciably heavier design. Also, the relatively long focal length and the 
positioning of the objective and observer’s cabin at opposite ends would call 
for a heavier connecting structure. All this points to a fairly easy choice, 
and further examination indicates the Cassegrain layout of reflector, since 
it offers both compactness and lightness. 


Optical Design 
The requirements of the main mirror may be set out as follows: 
(1) A high resolving power, especially if the observatory is to be used 
for visual observations of the Moon and planets. 

(2) Both photographic and very high power visual observations will 
require good light gathering capacity. 

(3) It must be readily transportable in a ferry rocket with accelerations 
of up to 12g. 

(4) Susceptibility to meteor damage must be as low as possible. 

(5) Whilst the useful life of an aluminized surface should be much greater 
in space than within our atmosphere, some measure of corrosion will 
still take place in the coating itself. It will therefore be advantageous 
to have a mirror construction that will permit re-aluminization to be 
carried out within the comparatively small space of an artificial satellite. 


An examination of the above requirements indicates a segmented mirror. 
By this I do not mean a tessellated mirror,? which has relatively poor resolving 
power, but rather a mirror having six to eight components all figured simul- 
taneously. Structural unity is afforded by a light alloy lattice, to which each 
glass segment is mounted in a manner permitting adjustment. For the purposes 
of grinding and figuring, the narrow spaces between segments are filled with 
glass sections and a plug is fitted into the central hole. The mirror is then 
worked and polished in the usual manner, and the filler sections subsequently 
removed. By supporting the mirror upon a network of pressure-balanced 
hydraulic jacks it could be safely transported in a ferry rocket. Should the 
question of size or floor strength prevent this from being done, it may be 
stripped down into its component parts and re-assembled in orbit. This would, 
however, require the complete readjustment of the mirror segments, which 
although theoretically possible would be difficult and very tedious. Never- 
theless, such adjustments will be necessary at infrequent intervals whenever 
the mirror is aluminized. Thickness of the segments must be sufficient to 
ensure that they retain their figure with passage of time, and to prevent 
permanent strain taking place during the high accelerations incurred in the 
ferry rocket. The usual problem of mirror distortion under its own weight 
at varying angles of elevation does not occur. 

Impact damage from a meteor would be isolated to the segment hit, and 
that could be removed or blanked off. It is unfortunate that it could not be 
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replaced, and the resolution of the mirror will certainly be impaired, but 
that is better than having no mirror. Moderate damage or distortion to the 
light alloy lattice may be overcome by readjustment of the remaining sections. 

Such a telescope would permit magnifications very much higher than any 
terrestrial telescope can approach simply because of the absence of any un- 
steady atmosphere. An image of the Moon at the prime focus of a 250 cm. 
(100 in.) mirror would be about 25 cm. (10 in.) diameter, and Mars at a 
favourable conjunction would present a disc 0-375 cm. (0-15 in.) in diameter. 
Photographs will record details which on Earth can only be observed visually. 


Structural Design 

In considering the structure of the telescope, it is of paramount importance 
that the main and secondary mirrors, together with the photographic plate, 
are rigidly connected to one another. They must also be shielded from small 
meteors and their connecting structure insulated against the effects of the 
Sun’s heat. 

The foundation of the observatory is the mirror lattice. Mounted from 
its front periphery and bore are three or four braced bipod cantilevers carrying 
the secondary mirror. From complimentary points on the rear face, thus 
obviating bending strain in the lattice itself, are attached the observer’s cabin 
and outer shell. The mounting for the photographic plate holder and its 
attendant instrumentation is fitted directly on to an extension of the lattice. 
By using this manner of construction, a number of difficulties can be 
circumvented. 

The observer’s cabin and outer shell prevent any stray light from entering 
the telescope and they also shield the major components from the Sun’s direct 
heat. The outer surfaces will probably be silver or aluminium in colour, but 
the inner surface of the shell will necessarily be black to prevent internal 
light reflection. Unfortunately, this means that the shell can more easily 
radiate heat across to the bipod cantilevers. This could be overcome by using 
a dull black lining of low heat transfer coefficient. A very light expandite 
product such as Onozote would serve admirably for this purpose. The main 
mirror, and the secondary mirror cantilevers, being independent of the outer 
shell, will be quite unaffected by the thermal distortions taking place in the 
latter. Similarly, a lack of stiffness in the shell and the cabin is permissible, 
as their inertia deflections will not impair the efficiency of the telescope. 

The mirrors are protected from all small meteors except those proceeding 
directly down the line of sight. These are unavoidable, but when the telescope 
is not in use it may be directed at the Earth, thus allowing that body to act 
as a shield. Generally speaking, it will be possible to operate the observatory 
right round the clock and the provision of mirror shutters is therefore hardly 
worthwhile. 

Basing calculations on the work of N. H. Langton,’ the shell and the 
bumper screen of the cabin would be constructed from chrome steel about 
1 mm. (0-040 in.) in thickness. In the case of the cabin, the bumper screen 
will cut down the temperature stresses in the aluminium pressure sphere and 
thus reduce the chances of fatigue failure. 
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For a telescope with a mirror of 250 cm. (100 in.) diameter, the total length 
of the observatory would be about 12 metres (40 feet). A project design reveals 
the following weight breakdown: 








Mirror, including lattice. . Ae ws .. 1,550 kg. (3,400 Ib.) 
Secondary mirror and supporting structure .. 270 kg. ( 600 Ib.) 
Outer shell ee ics - - -_ 770 kg. (1,700 Ib.) 
Observatory cabin - i - 540 kg. (1,200 Ib.) 
Air and power supplies, coma aquipasent a 450 kg. (1,000 Ib.) 
Astronomical equipment ry es a 230 kg. ( 500 Ib.) 

Total weight ag “ 4 .. 8,810 kg. (8,400 Ib.) 
Polar moment of inertia a 3-66 x 10’ kg.cm.? (40 ton ft.?) 
Lateral moment of inertia af 3-0 x 108 kg.cm.? (330 ton ft.*) 

Control 


Provision for control in three planes is required. For the purposes of 
telescope operation it is only necessary to provide turning couples, but the 
need to move the observatory relative to the spacestation calls for at least 
one mode of lateral acceleration. 

Three pairs of directionally opposed jets mounted on the control sphere 
are technically sufficient for obtaining initial orientation, but it would mean 
that one would occasionally have to turn the telescope through almost a 
complete circle. At the cost of a little increase in weight, it would be preferable 
to provide back-to-back jets. By temporarily disengaging the coupling 
system, they may then be used to provide translational thrust. Whether 
these jets are air, steam or gas does not immediately concern us. If each 
has a thrust of 0-45 kg. (1 lb.), giving a couple of 135 kg.cm. (10 Ib. ft.), the 
precessional acceleration is 4-4 x 10-* rad./sec.? 

To turn through a right angle, using jet power economically, one would 
maintain thrust for 10 seconds, rotate freely for 6 minutes, followed by a 
retardation period of 10 seconds. If we assume that the quickest manipulation 
possible is 0-25 second thrust immediately followed by 0-25 second reverse 
thrust, the angle turned through is nearly 6 seconds of arc. This is quite 
inadequate for the fine control necessary for photographic work and a secondary 
system is therefore required. 

Obviously electrical power will be available, probably in the form of batteries 
recharged at the spacestation. We can therefore utilize miniature electric 
motors equipped with small flywheels, in the following manner. Three motors 
are set up with their respective axes in each of the three planes. To steady 
the observatory after initial positioning, the motors slowly accelerate the 
flywheels until such time as the motion of the telescope is, in effect, absorbed 
by them. Subsequent corrections are made by successively accelerating and 
braking the flywheels. 
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To expect the observer to check three planes of motion and carry out 
the other functions required of him is probably asking too much. It is therefore 
likely that the two lateral planes at least will be switched over to photoelectric 
cell control. When a small meteorite hits the observatory, it would take a 
few moments for the flywheels to absorb the motion produced and for the 
telescope to be brought back on to its original line of sight. A relay switch 
will therefore have to close the shutter as soon as any one cell loses its guide 
star. In the event of appreciable movement having been incurred the photo- 
electric cell could be confused by the other star images. It is therefore better 
that the observer manually controls the re-opening of the shutter rather 
than allow it to be automatic upon the cells regaining apparent control. An 
impact leading to air leakage from the cabin would produce a continuous 
though decreasing torque. This would be beyond the capacity of the motors 
to correct, as a limitation is imposed by the maximum permissible speed of 
the flywheels. 


Pressurization 

The volume of the cabin is about 7-6 m.’ (270 ft.3), that of the airlock 
1-3 m.® (45 ft.’), and the internal pressure will be of the order of 0-7 kg./cm.? 
(10 lb./in.?)._ Such an observatory will be an expensive undertaking and we 
may assume that it will be kept in continual use by half a dozen or so 
astronomers. There is, therefore, little point in providing air reservoirs into 
which the air may be pumped when the cabin is unoccupied. Such a system 
may have reduced the leakage rate and lessened the chance of air loss by 
meteor impact. 

In order to keep the structural design simple and to facilitate maintenance 
of the mechanical equipment contained therein, the space beneath the cabin 
floor also would be pressurized. To ensure that the observer is given ample 
time to take the necessary corrective action, a pressure-sensitive warning 
system would be triggered off by any sudden increase in leakage rate. Should 
the pressure fall below a certain danger level, the reserve air supply will 
automatically discharge into the cabin and maintain the pressure. 


Operation 

As in the case of terrestrial observatories, the work of the telescope will 
be planned months in advance. Whilst only that hemisphere directly opposite 
to the Sun can be photographed at any one time, at certain periods it will be 
possible to photograph large areas of the star sphere with almost unlimited 
exposures. Although there will be no local atmosphere to diffuse moonlight, 
both the Earth and Moon will be troublesome if their light penetrates down 
the shell and is reflected thence on to the plate. When the Earth passes 
through known meteor streams the observatory would be placed downstream 
of the artificial satellite. With these considerations in mind, the possible 
duration of observation for each sector of the star sphere may be plotted 
over a suitable time base, say one year. With the aid of these charts and 
knowing the work to be accomplished, a suitable programme can be formulated. 
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When each astronomer goes across to the observatory he will take with 
him the plates required for the work he is to carry out. In the case of star 
fields or where any other large area is to be photographed, cut film some 
25 cm. (10 in.) square will be employed. These must be pressed to the correct 
shape in the focal plane, which on earth is normally done by vacuum from 
behind. This is obviously impracticable in the satellite telescope and to pro- 
vide positive pressure from in front is unattractive on a number of counts. 
The simplest solution appears to be to make use of an adhesive backing and 
press the film on to a formed block. The film can be peeled off prior to, or 
during the developing process. Low temperatures increase the sensitivity of 
emulsion (at the expense of contrast) but a few degrees above absolute zero 
is likely to be too low. In this event, the temperature of the film may be 
maintained by utilizing an electrically heated block. There being no air in 
the telescope, the plage of circulating air currents will not be realized. 

Upon entering the cabin through the airlock, the observer will check his 
air, fuel and electrical power supplies. Having divested himself partially or 
entirely of his spacesuit he must turn his attention to aligning the telescope 
on to his first objective. This is no simple matter as he has no accurately 
indexed equatorial mounting to work from. Probably using a wide-angle, 
low-power telescope he will determine its existing position (or check the 
figures given him by the immediately preceding astronomer). It will be 
convenient to retain the system of hour angle and declination and he will 
therefore first rotate the telescope so that the planes of control are laid closely 
along these reference lines. A quick calculation will then determine the 
durations of jet thrust and free swing in the two planes necessary to give 
the telescope approximately the correct aim. When there is a recognizable 
star or planet to centre on, this would now be used for the purpose. In many 
cases this will not be so and a sequence such as that indicated below will 
be required. 

Two suitable stars will have been chosen whose positions are accurately 
known, and their mutually subtended angle to the lines of declination calcu- 
lated. Centring on to one of these, the other is used to accurately align the 
polar planes of the telescope on to the hour and declination circles. By now 
using the motors, which are suitably calibrated for the purpose, the telescope 
is slowly precessed through the relatively perpendicular angles required. 

Two or three guide stars are picked out on the edge of the field and the 
photoelectric cells adjusted. The observer tightens his seat strap, opens the 
airlock in front of the plate, checks that the telescope is steady and opens 
the shutter. His cramped vigil will be similar to that of the astronomers 
down on Earth, but he will not have to endure the freezing night air, nor 
keep an eye open for clouds. 


REFERENCES 
(1) C. A. Cross, J.B.J.S., 1955, 14, 137. 
2) G. Horn D’Arturo, J. Brit. astr. Ass., 1953, 63, 71. 
(3) N. H. Langton, /.B.J.S., 1954, 13, 283. 
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PERTURBATION OF ELLIPTIC ORBITS 
BY ATMOSPHERIC CONTACT* 


By T. R. F. NONWEILER,f B.Sc., FELLow 


We here suppose a vehicle would be describing an elliptic orbit with focus 
at the centre of a spherical Earth, were it not that, at its perigee, it makes 
contact with the tenuous outer limits of the atmosphere (assumed at rest relative 
to the Earth’s centre). We shall suppose that the disturbance to the motion 
due to air resistance is sufficiently small that it may be treated as a perturbation, 
and we calculate the approximate change in the features of the instantaneous 
orbit. The approximation lies in the neglect of the effect of the perturbations 
during any one orbital circuit on the magnitude of the air resistance itself,! and 
certain assumptions about air density variation with height. 


Preliminary Analysis 

The vehicle drag D is supposed to be given by the law: 

D = 4pV°C,S (C, = constant) .. - ‘a (1) 
where V is the air speed, which in this context is also the speed relative to 
the Earth’s centre, S is some reference area, and C, the drag coefficient based 
on this area. The air density (p) at altitude (4) just above the perigee is 
assumed to be adequately described by an exponential variation with height: 
ss se = — bs (H, = constant) ae wh (2) 

pdah Hy “ 

This is strictly true only of an isothermal atmosphere of constant composition, 
in a uniform gravitational field, but in view of the rapid decline of density 
with height it is usually accurate enough to describe density variation over 
the limited range of height in an orbit where drag is most effective. 

The equation of an elliptic orbit in terms of polar co-ordinates (7,0) centred 
at the Earth, and with @ = 0 corresponding to the apse line joining the Earth’s 
centre with the perigee, is: 

r 1 — é 3 
a 1+ecos0 (3) 
where a is the major axis and e the eccentricity of the orbit. 

The speed in the orbit is 


2 1 + 2e cos @ + & 
y—4 = 1) =! v0 «tne Pos —_— 
avr a 1— é 


where (j/r?) is the value of the Earth’s gravitational acceleration. The distance 
moved by the vehicle in traversing an infinitesimal arc d@ is 


1 /dr\?|* a (1 — e*) 
= + = = ——,, (1+: s@ + e)t 
ds = rd6|1 b ¢ | (I + e cos 6)? (1 + 2e cos 6 + e*)*dé 


* Manuscript received 16 November, 1957, and in revised form 12 March, 1958. 
+ Department of Aeronautical Engineering, The Queen’s University of Belfast, Belfast- 
Northern Ireland. 
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dé a0 be ' - © 
whence: = § = Ered (l+ecos@)? .. rv (5) 


dt ds — 1 — é) 
The Perturbation from each Orbital Circuit 


In assessing the magnitude of the drag we ignore the change in the orbital 
characteristics during any one circuit, and take these as (say) corresponding 
to the instantaneous orbit at the perigee, distinguishing them by subscript zero. 


Thus, from (1), (2), (3) and (4): 


(1 + 2¢e, cos 6 + e?) ae the agg (1 — eg) sin? 4 0 
(1 + é)? ™ H, (1 + e, cos @) "i 





D=D, (6) 





where D, is the drag at the perigee. The change in velocity dV due to the 

action of drag in describing an infinitesimal arc d@ of the perturbed orbit is: 
ay ee 

~  m dO 

where m is the vehicle mass. Using here again the value of d@/dt appropriate 


- 


to the unperturbed orbit, we have from (5) and (6): 
dV ___ Df ag° (1 — &)*]# , éo? sin? 6 h 
d6-— m| pw (+ &) ~ (1 + eg) cos 6)? | 


= Whe Qagey (1L—“eg) sin? $0 7 
X xP] — Tr + cd) -» (7) 





From the results of the classical theory of instantaneous orbits,’ we note 
that a change in the speed dV is associated with corresponding changes da 
and de in major axis and eccentricity, and a rotation d@ (in the direction of 
the orbit) of the apse line, where: 








be re eV ad V 7% , 
Var" ind" Seta 


Once more replacing a and e by a, and e, where they occur in expressions 
multiplied by the drag (on the grounds of consistency), it follows that the 
changes Aa, Ae, and Aa, in a, e and @ per circuit, are given by: 
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The last relation, showing that the apse line does not rotate, follows from 
symmetry ; this reveals the deficiency of the approximation used, as, of course, 
the action of drag will not be entirely symmetric—the vehicle progressively 
losing energy whilst describing the circuit. 

To evaluate the integrals we make use of the fact that: 

H, Hy 
——- Eo me GCe |) .. as - 9 
Ay (1 — eq) R+Ny die 9) 
where R is the Earth’s radius, and hy the height at perigee. For we see 
from (7) that the change in speed at perigee per circuit is: 











2Q2r ” 
D s — £,)* |* 2 
Sicaives 2] (1 — ¢9) | | Fee Ta, Sint 40) d d 
m| p (l + é) l 
5 P .. (10) 
here fled.) = | Ue] A f =e 
where / (e, A, &) = (1 — ab? exp | (1 — Ab) J 


and for « 0, with 0 <A <1 and for any € in (0,1): 


SF (e,A, §) exp (— *) + O(e) bi . 


Neglecting the term of order e, it follows by using (11) in (10) and integrating, 


that: 
eae a) is(®) 
Vo 1 +e \ mg, 7° \2e ot - oe 
{ 


where jn (&) = exp (—&) J, (6&)/a" 





and similarly, from (8), we can calculate that: 


Aa, ) 
2--@s@) - - a 
Ae = — To Ae ee (st +in(3*)| . 6 


where A, = 2¢9/(1 + ¢) and ¢€, is given by (9). From (7) and (8) we may like- 
wise show the change in perigee height per circuit as: 


a (2)[a(8)-DIGED~ a 
a -iFs G o(3 h 2€, 1+ @/° 19) 


Thus the perigee approaches the Earth more and more rapidly at each circuit. 
The expressions J, are Bessel functions of the first kind, and their 
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behaviour for large values of Ag/€, (¢.e., for orbits in which the difference between 
heights at apogee and perigee is large compared with H,) is given by: 


a) = a0(-2)0(B)~CC = re 
(Babeo(-2)a(B)()0-22-BE--) 


whilst for orbits of such small eccentricity that eg < €, (¢.e., the difference in 
heights between apogee and perigee is small compared with H,) then the 
following expansions are useful: 


do Ly SAP 
Yr -)= “Fa hwy 


> 
mY (Ie, SM) 
2€5 2 | Wee A | J 


Perturbations over a Number of Orbital Circuits 

We now treat the perturbation from each circuit as infinitesimal, which 
from (13) and (14) either implies the same assumption about the magnitude 
of the parameter ¢, defined in (9) or the parameter (D,/mg,). Further, we 
shall suppose that €, is the same for each circuit. This may, in fact, be 
objectionable in practice, particularly as the height at perigee decreases pro- 
gressively. However, it is not difficult to generalise the results we shall give 
so that a different value of ¢ is taken for consecutive groups of circuits, over 
which the height of ry perigee is within stipulated limits. 

From wi and (14) we have: 


“C3 er *) | 20e Jo (Ag/2€o) — 1 (Ao/2€0) Ja. 
R+ Tk Cojo (Ag/2€o) + Ja (Ao/2€o) 
Summing each side over a number of circuits, by treating A as an operator 


d/dN and integrating with respect to N, = number of circuits, we find for 
small (constant) e that, denoting by suffix “7” the initial conditions at perigee, 


é;) ie Ah, (A A 
(<)m ioe = fies pei (Ao/2e) is oe 


the error terms omitted from the right-hand side being of order ¢. From (2) 
and (9) we see that: 


= [a=] -O* GR ee 


the approximation here being valid since the range of perigee heights over 
which air resistance has any appreciable effect on the orbit is small compared 
with R. Quoted in this form, the left-hand side of (18) would be essentially 
independent of the choice of H, (or of €), and from (16) and (17) so also is the 
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right-hand side, unless the eccentricity changes from a finite to a small value 
(of order €) which is the only condition leading to the ratio of the modified 
Bessel functions being sensitive to e. 
A similar argument to that used in devising (18) shows that, using (12) 
and (14): 
V2 14+ 64 


VA i1+e,' 





amen? | oe 


which (to this accuracy) is in agreement with the expression for the perigee 

speed of unperturbed orbits (4), as the height at perigee a, (1 — eg) changes 

only by a fraction proportional to «. From (1), (19) and (20) it follows that: 
Dy — Dy & (1 + 0)? J, (Ay/2€) 


mM§o mg; € (1 + e)* J; (Ag/2€) 


with error terms omitted of order « compared with those preserved. Thus 


n (14): 
a a) (i + te , odo (Ao/2¢) (Ao/2€) (A,/2¢) 
= _ mg;/ € (1 + @;) Y FulAg/ 2) Ja (Ail 


On integrating as before, if Nj» is the number of circuits involved in changing 
from the initial perigee eccentricity e; to é), then: 


r (1 + e,) mg; ey 
Nig = ane, nt 4me, j, (A;/26) 55 (7 >) Fes + + In ( 1+ | any 


where we have once again ignored terms of order € compared with unity. 
More simply, using (1), (4) and (20) we have alternatively: 


Nig = N(e;) — N(eo) 

(m/C pS) . e 
Imp (R+H|™O+9-TT, 
and where, from (18) and (19), the density at any perigee, including the initial 
one, is given by pp where 











(21) 
where N(e) = 


we 1 (A/2€)p9 = p*, a constant ae -- (22) 


Then equations (21) and fe provide all the information necessary to compute 
the orbital characteristics e, and a, at any circuit. In (21) we note that for 
small values of e 


é os a 
in(i +9 —i7, 73° as 4? 


For some purposes it is desirable to relate the changes to the passage of 
time, and in the present notation, from equation (5), the periodic time per 
circuit is: ten! 
=T1Ag 
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on ignoring the changes in the instantaneous orbit due to air resistance over one 
circuit. The total time lapse ¢,, in describing Nj» circuits can be found by 
a similar argument to that used in obtaining Nj», and the result is: 


tio * t(e;) ——- t(éq) 


tt(e 
tere —te)_ —_ me) : 
where (R fi hy? 4p* (R j h) [7% — 70) | 
3+e 1 3—e 
: Tl. Pes te Te bas ee eB 2 
and 7(e) d+ed—ert~ cosh G i" :) 2 . 


which like (21) is an expression correct to the first order in €. For small 
values of e, we note that: 


re 7  . 
T (e) constant + 2e? < — 3¢ +... 


The equations of this section become indeterminate if the initial orbit is 
circular (so that e = 0). From (14) there is then, of course, no change in 
eccentricity arising from the perturbation due to air resistance; but, the radius 
of the circular orbit reduces and, in fact, from equations (12) to (15), putting 


= 0 we have: 
2\V Aa Ah dw >) 
Va J my 


The neglect of changes in h compared with R in computing the drag is then 
equivalent to the neglect of changes of speed in the circling orbit. Either 
by direct integration of the above expressions, or by proceeding to the limit 
e — 0 in (21) and (24), it can be shown that the number of circular orbits is: 


re (m/C,S) 5 1 
. i 2a (R+ Pi —5) 


and the time involved is: 


tig = (m/C,S) (A=) (- o ~)s 


Summary of Analysis 


The perturbation on a single circuit of an elliptic orbit has been evaluated 
neglecting the effect of this perturbation on the magnitude of the air resistance, 
treating the atmospheric density as varying according to the exponential 
law (2) at altitudes just above the perigee. The essential results are contained 
in equations (13) and (14) giving the changes in major axis and eccentricity 
on the circuit. The notation is summarised at the end of this paper; the 
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modified Bessel functions occurring are tabulated, for example, in Ref. 2 and 
their variation indicated in Fig. 1 and 2 and by equations (16) and (17). 
We next assume that the air density variation can be approximated by 
a number of piecewise exponential variations, so that in (2), the value of H, 
has a different value over different altitude ranges. Then (provided these 
height ranges are small compared with R and the drag coefficient at perigee 
conditions does not vary) we can work out the aggregate perturbation for a 
number of circuits such that the initial and final perigee heights fall at, or 
within, the limits of such an altitude range. The aggregate perturbation over 
a larger range of altitude can of course be evaluated by repeated application 

















0-8 ad | 
. ae 





—— — = 
























































0 | 
w | | 
yer | N_| 
0-6 | | | h 
| | | 
0-5 
\* 
at | | — 
at | | 
“ | | ae ig | | 
= | ACD NX | 





0 2 + T \ 
| eal 












































0-1 
iT | | | 


0-1 0-2 0-4 0-6 0-8 | 2 at es 
g 


Fic. 1. Variation of modified Bessel functions used in the analysis 














PERTURBATION OF ELLIPTIC ORBITS BY ATMOSPHERIC CONTACT 375 


0-9} —J [ —+—_+— vain 
/ 























j(§)Hi(8) 


0-6 \ —+—_}— 














a 





0-4 | " 





0-3/-——_|—+- 


; | 7c 


0-1 0-2 0-4 06 0-8 | 2 4 6 8 10 


j(€)-i,(€) 
~*~ 
rd 

































































Fic. 2. Variation of modified Bessel functions used in the analysis. 


of the results, piecewise. The essential results are the relation between 
eccentricity and the number of circuits, given by equation (21), and that between 
the height at perigee (as determined by the local air density pg) and eccentricity 
in (22). The time scale of the process is given by (24). 

In all examples a knowledge of the drag coefficient is required: usually 
for bluff bodies moving at the speeds involved in elliptic orbits a coefficient 
of 2 on frontal area is expected. At the time of writing no reliable information 
is available on the density at great heights, but the data of Table I derived 
from Ref. 3 and 4 give a tentative estimate. 
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TABLE I.—Density Variation with Height 


























I 
Density Density | 
relative to relative to 

Height that at Hy Height that at i, ] 

(miles) Sea-Level (miles) | €, x 10°|} (miles) Sea Level (miles) | «, x 10° 
0 l 5:3 1-35 120 *2-61 x 10-* 149 | 3-65 
10 1-33 x 10-4 4-0 1-01 130 1-40 x 10-1 17-1 4-2 
20 1-10 x 10°? 4-3 1-11 140 8-0 x 10-4 19-3 4-7 
30 1-18 x 10-3 5-1 1-25 150 49 x 10-0 21-5 5-2 
40 1-68 x 10-4 5-0 1-25 160 33 x 10 23-7 5-75 
50 1-66 x 10-5 3-8 0-95 170 21 x 10-4 259 | 63 
60 1-23 x 10-* 4-3 1-07 180 1-46 x 10-4 28-1 | 68 
70 1-19 x 10-7 5-3 1-31 190 1:03 x 10-" 30-3 | 7:3 
80 1-79 x 10-8 7-1 1-76 200 76 x 10-8 343 | 82 
90 4-21 x 10-° 8-8 2-18 210 58 x 10-8 38-3 | 92 

100 1-28 x 10-* 10-5 2-56 220 45 x 10-® 42-3 | 10-1 
110 5°39 x 10-° | 12-7 | 3-10 230 36 x 10-8 46-3 11-0 











Note.—Data for 20-100 miles height from Ref. 3; interpolated with data relating to 
heights above 150 miles given in Ref. 4. Sea Level density = 0-0765 Ib. /ft.*. 

* Provisional data for Satellite « (1957) suggests that this figure is less than half the 
true value at this height. 


Application of Results 


There are three problems which are of interest in the application of the 
results we have derived: these arise from the three separate physical processes 
involved, namely, the orbital geometry, the nature of the atmosphere, and 
the aerodynamic forces on the vehicle (relative to its weight). If any two of 
these are completely specified, then something about the third may be derived. 

We have so far stated the problem as if the interest lay in the determination 
of the orbital geometry, but it might lie in the determination of the density 
variation with height. We see from (13), for example, that the change in 
major axis—which, of course, may be related to the change in periodic time 
(49) by equation (23) as: 

At, 3 Aa, 


25 
ty 2 ay aaa 

is determined not only by the drag of perigee (Dy), and so by the air density 
there, but also by the exponential rate of density decrease (defined in equation 
(2) by the height H,). An exception arises if the orbit is so nearly circular 
that the density on a complete circuit is more or less the same. Using (13) 
and (16) we see, in fact, that, for differences in heights at perigee and apogee 
which are large compared with Ho, then 

Aa 

—* Dect & polo? 

4 
In other words, so far as the atmosphere is concerned, not only the density 
at perigee (p,) but also its variation at greater heights is important. Thus 
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two items of information are necessary to specify the atmospheric density. 
Suppose, for instance, we also knew the rate of approach of the perigee. 
Now from (13) and (15): 


Ah, _ (1-% hr (Agl2¢) 
May “G ey [1-2 Jo (Ao | 


oy re 
if Zaye, > Hy 





a5 


the approximation here arising from using (9) and (16). Thus we have a 
direct determination of Hy. We see, in fact, that the apogee approaches 
the Earth at a rate (2a9¢,/H,) times faster than the perigee (the distance of 
heights between perigee and apogee being twice the focal distance age). 


TABLE II.—Typical Orbit of Small Eccentricity 











| Height | Height | Aa Ah, Ato /te 
at at (feet (feet (sec- Cir- Time 
Peri- | Apo- | Period per per onds cuits | from 
Eccen-| gee gee (min- cir- cir- per Ae (per from | launch 
tricity | (miles) (ates utes) cuit) cuit) day) circuit) launch (day s) 
0-05 125 555 95-5 223 83, 1:28 |9 x 10-* o 0 
0-04 123 463 94-0 270 124; 1-55 | 1-1 x 10-*| 1,009 66 
0-03 120 372 92-4 350 21 2-03 | 1-4 x 10-5 | 1,805 118 
0-02 117 283 90-9 510 46 3-03 | 2-1 x 10-* | 2,360 153 
0-01 111 193 89-7 | 1,070 186 6-36 |4 x 10-*| 2,713 175 
0-005 104 145 88-45 | 2,590 880 | 15-5 7-9 x 10-* | 2,807 181 
0-002 924 109 87-7 | 12,600 | 7,500 | 75-8 2-3 x 10-* | 2,834 183 
0-001 85 93 87-3 | 43,600 32,900 5-1 x 10-* | 2,837 183 
0-0005 79 83 87-0 On the last | circuit 2,838 183 




















To illustrate the direct application of the results to the determination of 
orbital geometry, an example has been concocted in which a vehicle, whose 
air resistance can be computed assuming (m/C,S) = 25 lb./ft.?, is imagined to 
be launched into an orbit of eccentricity 0-05 at the perigee of 125 miles altitude. 
It gradually slows down in an atmosphere whose characteristics are supposed 
to be identical with those stated in Table I. The details of the orbit are then 
given in Table II, a different constant value of the parameter e being assumed 
over each of the intervals between successive conditions enumerated therein. 
It will be observed that only over the last few circuits does the air resistance 
produce any considerable perturbation to the orbit, and so the results are 
valid except possibly to the dying moments of the lifetime of the vehicle, 
which it will be seen is about 6 months. The total lifetime may, with good 
approximation, be taken from equation (24) as simply ¢(e,), which works out 
as equal to 179 days, as compared with the more accurate figure of 183 days 
derived from the step-by-step procedure (using different values of «). Fig. 3* 


* Not reproduced (see final paragraph). 





Te ayes =F." 


“wn: 


A ~~ ON - e 





378 T. R. F. NONWEILER 








is obtained using this approximation, together with the following formule, 
thereby derived from (13), (23) and (25): 


te) 3(1 +e)? (1 —e,)"* 75(A,/2€) 

= wa . (Tle) — TO] .. 26 

t 8e,( — At;/t,) 9,(A;/2€) Me) — 70) (9) 

where (—A?,/t;) is the initial time rate of decay of orbital period. In the usual 

condition in which e;>e, the series expansion of (26) in ascending powers of e; 
and (e/A,) yields: 


t(e;)At; _ 3 ati 7 ae 
2 =e + gate) +yt--); 


For near parabolic orbits, a good approximation to the lifetime is given by 
¢(1), where: 











t(1)At, 

Of course, such orbits might be perturbed by bodies other than the Earth. 
Otherwise, if ¢; is of order e, then from (26): 

004, _ 3 eel 

t2 45 (ele) 





= — 3[1 + O(6)} 





(1 + O (e)). 


In particular, for an initial orbit which is circular (e; = 0), its lifetime is 
(0), where: 








(0) 3 € € mg; 
9 


a well-known and simple result, though one of little utility as orbits are unlikely 
to be so closely circular that its use is justifiable. 

According to the analysis—and in particular equation (22)—the orbit does 
not become circular before ground impact: in other words, the vehicle eventually 
descends to the ground because of the approach of the perigee rather than a 
decrease in eccentricity. The last circuit has still a slight eccentricity, though 
its perigee is below ground level. Thus the final descent is virtually at a con- 
stant speed, rather above the circling speed, the air resistance determining the 
angle of descent. This is perfectly compatible with the results of King—Hele’s 
analysis of descents from circling orbits,® and, indeed, it will be seen in Table II 
that the penultimate circuit is almost circular. Of course, it is to be expected 
that, in actual fact, the speed will drop below the circling speed in the final 
circuit, because neither the present analysis, nor that of Ref. 5, are applicable 
where the drag retardation is so large as it would be in this ultimate penetration 
of the lower atmosphere. 

The final descent is quite precipitate: even at the perigee of the penultimate 
circuit the equivalent air speed is only about 3 ft./sec. This fact is still generally 
valid however large or small the drag “‘loading”’ (m/C,S) happens to be. Thus 
the aerodynamic problems of atmospheric re-entry are not appreciably alleviated 
by the protracted period of gradual deceleration of the elliptic orbit. 
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In conclusion, the author wishes to thank Dr. D. C. M. Leslie of the Royal 
Aircraft Establishment, Farnborough, for indicating an error in the treatment, 
which has been corrected in the text presented here. However the correction 
was made too late to allow modification to Fig. 3, which has therefore been held 
over until a later issue of the Journal, where it is hoped some further deduc- 
tions from the present analysis will be given (including some results deduced 
from the first satellite launchings). 


LIST OF SYMBOLS 


Cp drag coefficient of vehicle based on frontal area. 


D drag. 
H, inverse logarithmic rate of density change with height (defined by equation (2)). 
N number of circuits. 


Nj number of circuits from initial to present perigee. 
R Earth’s radius. 


S frontal area of vehicle. 

T(e) defined in equation (24). 

V speed of vehicle. 

a instantaneous semi-major axis of orbit. 

e instantaneous eccentricity of orbit. 

g gravitational acceleration. 

h height from Earth’s surface. 

jo(&) modified Bessel functions of — defined in (12); Fig. 1 and 2, and equations 
dx(8) (16) and (17) aid computation. 

m mass of vehicle. 


distance of vehicle from Earth’s centre. 
s distance measured along orbit. 
t time. 
tio time lapse from initial to present perigee. 
to periodic time. 
A operator denoting rate of change per circuit. 
€o = H,/(R + hy). 
6 angle subtended at Earth’s surface between vehicle in its orbit and perigee 
(true anomaly). 


A = 2e/(1 + e). 
i” constant of gravitation. 
@ longitude of nearer apse. 
p air density. 
p* constant density defined by equation (22). 
Suffix “‘0’’ denotes instantaneous conditions at present perigee. 
“= - - mm » initial perigee. 
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B.I1.S. NEWS 


B.LS. NEWS 


Space Medicine Symposium 


The Society is organizing an international symposium on space medicine 
and biology. It is to be held in London early in October, 1958. 

Future announcements will give the detailed programme and particulars of 
the arrangements for registration, obtaining hotel accommodation, etc., but 
members and others who contemplate attending the symposium can obtain 
further information from the Secretary. It is hoped to arrange the dates of 
the symposium so that participants will also find it convenient to attend the 
celebrations of the twenty-fifth anniversary of the founding of the Society. 

Although it does not receive the same publicity as the launching of artificial 
satellites, a considerable amount of work on space medicine is now being carried 
out in various countries, and the outcome of this work is of vital importance 
to the future development of astronautics. Especially noteworthy are the 
research programmes of the U.S. Air Force Department of Space Medicine and 
other American institutions, which include experiments on weightlessness, crew 
behaviour during interplanetary voyages, food and air supplies, and the effects 
of cosmic radiation (a particularly courageous investigation was the MANHIGH 
flight of Major David G. Simons—-see illustrations opposite). Members will 
already know about the Russian use of dogs in Sputnik II and upper-atmosphere 
rockets. In Britain there is the work of the Royal Air Force Institute of Aviation 
Medicine, and there are similar institutes in other European countries. 

It is hoped that speakers from several of these bodies will present papers at 
the symposium. The last date for the submission of papers is 30 June, 1958, 
but intending speakers are asked to write as soon as possible to Dr. A. E. Slater, 
Dell Farm, Whipsnade, Dunstable, Beds., advising him of their proposed 
subject. 


Election of Members 


The following elections were made at the Council Meeting of 4 January, 1958: 
Fellows. 
Joun RussELt Apams, 210, Princes Avenue, Palmers Green, London, N.13. 
KENNETH HENRY BLAck, 25, Vicarage Road, Chellaston, Derby. 
Joun Tuomas B Lake, B.S., 486, Orange Avenue, Merritt Island, Florida, U.S.A. 
ALAN BriGurty, “Electra,’’ High Road, Abridge, Romford, Essex. 
CARL-JOHAN CLEMEDSON, M.D., Vintervagen 35, Solna, Sweden. 
KEITH ALLAN CORKILL, M.Sc., 28, College Crescent, London, N.W.3. 
MICHAEL GEORGE HARTLEY, 208, Burton Road, Derby. 
FRANCIS J. HEYDEN, Ph.D., Georgetown University, Washington, 7, D.C., U.S.A. 
Ian Marcus Hossy, B.Sc., P.O. Box 64, Sabie, East Transvaal, South Africa. 
Puitip LEIGHTON, 2, Greenacre Road, Woolton, Liverpool. 
BENN DEVEREAUX Martin, B.S., JPL, 161-225, 4800, Oak Grove Drive, Pasadena, 
California, U.S.A. 
Joun DEREK Povan, 24, Copes Way, Chaddesden, Nr. Derby. 
Joun Rupkin, B.Sc., 417, Brereton Avenue, Cleethorpes, Lincolnshire. 
KEITH SHERWIN, Rose Bank, Kirk Ireton, Derbyshire. 
GEORGE CHARLES SZzEGO, Ph.D., Fallis Road, Loveland, Ohio, U.S.A. 
CHARLES TILGNER Jnr., B.S., Grumman Aircraft Engineering Corp., Bethpage, N.Y., 
U.S.A. 
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CLIFFORD WILFRID Tinson, F.R.A.S., “Three Winds,’’ Gloucester Road, Cheltenham. 
EUGENE EMILE VIELLE, R.A.F. Club, 128, Piccadilly, London, W.1. 
Joun WittiaMs, 15, Jefferson Drive, Brough, East Yorkshire. 


Senior Members. 


FRANK ALFRED Barr, B.A., 27, Harrow Road, Leicester. 

RONALD VINCENT BAUGHAN, “‘High Oak’’, 4, Hurst Green Road, Hurst Green, Oxted, 
Surrey. 

ANTHONY CHARLES CLARKE, 9, Wellesley Court, Fitzalan Road, Littlehampton, Sussex. 

ALAN MILNEs CrosTHWAITE, 4, Osborne Road, Cale Green, Stockport, Cheshire. 

GEORGE Brian Davis, 197, Bedford Hill, Balham, London, S.W.12. 

GILHAM ROBERT GrBBs, 37, Denziloe Avenue, Hillingdon Heath, Middlesex. 

Joun GorpDon, 36, Fairfield Avenue, Ruislip, Middlesex. 

MICHAEL JAMES HENDRIE, “‘Anchorage,’’ 25, The Crossways, Westcliff-on-Sea, Essex. 

ALFRED WILLIAM LuckKHuRsT, 469, Burgess Road, Swaythling, Southampton, Hampshire. 

HAROLD MartTIN, 20, Marwood Drive, Barnard Castle, Co. Durham. 

EpwarpD ALBAN TISDALE, 80, Dartmouth Avenue, Walsall, Staffordshire. 

Joun TaLsBot ULLMAN, 64, Rosebery Road, Epsom Downs, Surrey. 


Members. 


JoHN REGINALD Baker, c/o National Provincial Bank Ltd., 33, Kings Road, Chelsea, 
London, S.W.3. 

Davip Dean BEYER, M.S., 2611, Otis Lake Road, Anchorage, Alaska. 
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OBITUARY 


Christopher Rodney Armstrong 

C. R. Armstrong, an active member of the Society, died suddenly on 
5 November, 1957, at his home in London. 

After the war, during which he served in submarines, he studied at Man- 
chester University, graduating B.Sc. in 1948. He then came south to Guy's 
Hospital to study for the degrees of M.B., B.S. It was while he was there that he 
joined the Society in 1950. During the next seven years he took a great interest 
in the activities of the Society, helping in the preparation of the last membership 
list and lecturing to many outside bodies on its behalf, including the 1954 
Conversazione of London University, at Cambridge, and at Guy’s itself. 

During his medical studies he had switched to dental medicine, obtaining 
his L.D.S. in the June of last year. He returned to the last part of his medical 
course expecting to graduate in 1958. His interest in many aspects of natural 
science was insatiable (he was a Fellow of the Royal Astronomical Society and 
a Fellow of the Royal Meteorological Society). He conducted, with hospital 
approval, his own research on dental malformation, and assisted some of his 
tutors in their researches, which included investigations into the growth of 
crystals. 

There seems little doubt that had he lived the mass of information he was 
acquiring and his great enthusiasm would have enabled him to have made a 
valuable contribution in the field of space medicine. 


PERSONAL NOTES 


Mr. T. R. F. NoNWEILER (FELLOW) has left the College of Aeronautics, 
Cranfield, and has been appointed Senior Lecturer in the Department of 
Aeronautical Engineering, of the Queen’s University of Belfast. 

Dr. A. E. Roy (FELLow) has been appointed Lecturer in Astronomy at 
Glasgow University, with effect from April 1, 1958. Dr. Roy is Chairman of 
the Society’s Scottish Branch. 


P. R. F. 
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ABSTRACTS 
Edited by J. HuMpuHRIEs and J. FoLey 


A list of abbreviations of titles of journals was included with the 1954 
index and addenda have been published in subsequent issues of the Journal. 


ASTRONAUTICS 


(263) The “Clock Paradox” and space travel. E. M. McMillan. Science, 126, 
381—4 (Aug. 30, 1957). A résumé of the theory giving examples of the time dilatation 
effect. (3 refs.) 

264) The second artificial Earth satellite. Nature, 180, 931 (Nov. 9, 1957). 
Editorial comment. 

(265) Rocket of the first Earth satellite. Nature, 180, 944 (Nov. 9, 1957). Analysis 
of photographs taken at Dunsink Observatory, Ireland, gives orbital period of 95m. 53s. 
at Oct. 11-23d. decreasing uniformly by 4-13s. a day to 94m. 49s. at Oct. 26-78d. The 
corresponding decrease in semi-major axis is 3-4 km. a day. 

(266) Further radio observations of the first satellite. I. R. King, G. C. Mc- 
Vittie, G. W. Swenson and S. P. Wyatt, of University of Illinois. N. Lea of Marconi’s 
Wireless Telegraph Co. Ltd. Ministry of Supply. Nature, 180, 943-4 (Nov. 9, 1957). 
Three letters describing methods and observations. Univ. of Illinois gives period during 
interval Oct. 5-24 as: 

P = 96m. 09-8s. — 0-1293 s. n — 7-797s. x 10-5 n*? where n is the number of revolutions 
since an adopted zero, which was the southbound crossing of their parallel on Oct. 5d. 
13h. 30m. Ols. U.T. 

(267) Radar observations of the first Russian Earth satellite and carrier 
rocket. Staff of the Jodrell Bank Experimental Station, University of Manchester. 
Nature, 180, 941-2 (Nov. 9, 1957). Brief description of equipment used—no results. 
(3 refs.) 

(268) Observations on the orbit of the first Russian Earth satellite. Staff of 
the Royal Aircraft Establishment, Farnborough. Nature, 180, 937-41 (Nov. 9, 1957). 
Results of observations extending over most of the three weeks during which the satellite’s 
transmitter was functioning. Attention was focused on estimating the orbit constants 
although many data remain to be analysed. Estimated orbital data for noon G.M.T. on 
Oct. 15 are: 

Height going north at latitude 51-2° N 

Height going south at latitude 51-2° XN. 

Inclination of orbital plane to equator, assumed to be 65°. 

Perigee occurs going north at latitude .. +o ae 

Perigee distance from Earth’s centre 

Height at perigee over mean Earth 

Apogee distance from Earth’s centre 

Height at apogee over mean Earth 

Eccentricity es we 

Semi-major axis 

decreasing at .. - 

Harmonic mean radial distance - 

Period + ; 

decreasing at 


Rate of rotation of orbital plane (retrograde) 
increasing at 0-0029 deg. per day. 


Rate of rotation of major axis (theoretical) “ — 0-40 deg. per day. 

Track crosses latitude 51-71° N. at longitude 0-7° E. at time , 21hr. 08min. 56s. G.M.T. 
(2 refs.) 

(269) Flying tothe Moon. Engng., 184, 642-3 (Nov. 22,1957). Extensive abstracts 
of two Russian papers: ‘‘Problems in the dynamics of lunar flight’’ by V. A. Yegorov deals 
with computer calculations of free orbits, the effects of initial errors and the conditions 


130 + 6n. miles 
250 + 5n. miles. 


40-9° N 

3,561 + 7 n. miles. 
123 + 7 n. miles. 
3,922 + 8n. miles. 
484 + 8n. miles. 
0-048 + 0-002 
3,741 n. miles. 

1 n. mile per day. 
3,733 n. miles. 

95 m. 48-5 s. 

2-28s. per day. 
3-18 deg. per day. 
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required for capture by the Moon. ‘Some variation problems connected with the launching 
of artificial satellites of the Earth,’’ by D. E. Okhotsimskii and T. M. Eneev investigates 
the most economical flight path and thrust programme (printed in full, /.B.J.S., 1958, 16, 
263-294). 

(270) Satellite review. 1. The historical perspective. G. J. Whitrow, 
2. The political consequences. W. Elliot. 3. The strategic significance. A. 
Buzzard. Impulse, (3), 5-11 (Dec. 1957) A general non-technical review of the various 
aspects. 

(271) Life-time of an artificial satellite. J. A.Fejer. Nature, 180, 1413 (Dec. 21, 
1957). A simple method of obtaining the approximate life-time from initially measured 
values of the rate of change at apogee and the difference between the heights at apogee 
and perigee. 

(272) Estimating the life of a satellite. J. M. C. Scott. Nature, 180, 1467-8 
(Dec. 28, 1957). A simple method of obtaining the life, given the period and some know- 
ledge of the eccentricity—air density is not needed. 

(273) Scientific observations of the artificial earth satellites and their analysis. 
H. S. W. Massey and R. L. F. Boyd. Nature, 181, 78-80 (jan. 11, 1958). A report on the 
Royal Society meeting on the results of observations of the first two Russian artificial 
satellites. 

(274) A photographic observation of the satellite 1957 Beta leaving the Earth’s 
shadow. H. Neckel. Nature, 181, 257-8 (Jan. 25, 1958). Calculates co-ordinates 
from a photograph taken on Dec. 15, 1957 in the U.S. of the second Russian satellite. 


BIOLOGY AND MEDICINE 
(275) A technique for instrumenting sub-gravity flights. G. J. D. Schack and 
D. G. Simons. J. Aviation Medicine, 28, 576-82 (Dec. 1957). Apparatus for aiding the 
pilot in obtaining the maximum duration of the weightless state and to record the forces 
acting on the subject and aircraft. (4 refs.) 


PROJECTILES 
(276) Guided missiles and rockets. A bibliography, 1946-1956. Engineer 
School Library, Fort Belvoiy, 50 pp. (1956.) 
(277) Choosing a missile guidance and control system. R. W. Mayer. Elec- 
tronic Equipment, 5, 36—9 (July, 1957). Analysis of pertinent variables. 


RADIO AND ELECTRONICS 
(278) Observations in Australia of radio transmissions from the first artificial 
Earth satellite. G. H. Munro and R. B. White. Nature, 181, 104 (Jan. 11, 1958). 
Information on the fading of the 20 Mc/s signal. 


ROCKET MOTORS 


(279) Propellant vaporization as a criterion for rocket engine design: calcula- 
tions of chamber length to vaporize a single n-heptane drop. R. J. Priem. 
N.A.C.A. Tech. Note 3985, 41 pp. (July, 1957). Calculations based on droplet-evaporation 
theory show that for a given chamber length the percentage of fuel vaporized can be 
increased by decreasing the drop size and the initial drop velocity, or by increasing chamber 
pressure, final gas valocity and initial fuel temperature. 

(280) Chemical rocket propulsion in the Vanguard satellite launching vehicle. 
J. E. Burghardt. Chem. Engng. Progr., 53 (10), 86, 88, 90 (Oct., 1957). General description 
of General Electric X-405 first-stage motor using liquid oxygen and kerosine, and Aerojet- 
General AJ10—-37 second stage motor using white fuming nitric acid and unsymmetrical 
dimethylhydrazine. 

ROCKET PROPELLENTS 

(281) Processing rubber-base composite rocket propellant. C. F. Dougherty. 
Chem. Engng. Progr., 53, 489-92 (Oct., 1957). Describes preparation of materials, mixing, 
extrusion, trimming and final assembly in rocket case. Deals only with propellents 
containing ammonium nitrate as oxidant. 

(282) Anhydrous hydrogen peroxide as a propellant. R. Bloom and N. J. 
Brunsvold. Chem. Engng. Progr., 53, 541-7 (Nov., 1957). Properties are discussed and it 
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is compared with other propellants (inc. 90% peroxide) from a performance point of view. 
(9 refs.) 

(283) Possible fuels for Sputnik II. Engng., 184,.644 (Nov., 22, 1957). Suggests 
that hydrazine or pentaborane, with liquid oxygen, may have been used to launch Sputnik 
II. 


REVIEWS 


The Amateur Astronomer. By Patrick Moore. 84 x 5}in. 317 pp., with 61 figures, 

14 maps and 16 plates. 1957. London: Lutterworth Press. (25s.) 

Patrick Moore’s theme, in this fine new book, is to show the amateur astronomer, whose 
equipment is nothing more than a small telescope, just what he can do. After four intro- 
ductory chapters, chapters have been devoted to the Sun, the Moon, Occultations and 
Eclipses, Aurore and the Zodiacal Light, the Nearer Planets, the Outer Planets, and 
Comets and Meteors. A further seven chapters deal with stella astronomy in very clear 
terms. 

There is a brief chapter on Telescopes and Observatories, but the prospective buyer 
must not expect to find here adequate instructions for making a telescope; the omission is 
perfectly justifiable in view of the fact that there are numerous books which cover instru- 
mental problems. Instead, the reader is referred to a bibliography, forming the last of 
twenty-nine extremely useful appendices, which add considerably to the value of the book. 

Mr. Moore’s style is lucid, and he often introduces apt analogies. On occasions, the 
reviewer was amused by witty footnotes, but it must not be thought from this that the 
subjects are treated at all facetiously: the reading is light but, at the same time, it is 
crowded with useful facts. 

Among the appendices, one finds a very clear, simple map of the Moon, in two halves, 
each half facing a photograph of the same part, and there is an associated guide to many 
selected formations. There are fourteen clearly labelled star charts, including two key 
maps, which will enable the inexperienced observer to find his way about the sky very 
easily. Notes accompany each chart. 

Mr. Moore has considered as “‘completely false’’ Sidgwick’s statement that “‘if the 
aperture exceeds about 12 inches, the atmosphere will seldom allow the full aperture to be 
used’”’ (page 113). Certainly, the larger the telescope the more detail will it show on 
occasions when the seeing is good, and this seems to be Moore’s point. Sidgwick’s statement 
refers, however, to the amount of light which is collected, and, in his statement, the word 
“‘seldom”’ is important. For example, if the seeing is better than one second of arc, an 
aperture of 20 inches may perhaps be used without loss of light. A few other incorrect or 
misleading statements are to be found scattered throughout the book, but any book which 
deals with debatable problems should be read with an open mind. 

The word “‘parabola’”’ in Fig. 47 (page 137) should read “‘ellipse’’. Unimportant mis- 
prints were noticed on pages 126, 142, and 192. 

In addition to many line drawings, there are 16 full-plate photographs which add appeal 
to Mr. Moore’s well-compiled book. There is a good index, and, without doubt, no keen 
amateur astronomer should be without this very moderately-priced reference work. 

G. FIELDER. 


Sputnik into Space. By M. Vassiliev, under the supervision of Professor V. V. 
Dobronravov (originally published under the title IWYTEIMECTBUA B KOCMOC, 
translated by Mervyn Savill from the Italian ‘“‘Su Sputnik nei Cosmo’’. 8$ x 5} in. 
Pp. 147, with 19 illustrations. 1958. London: Souvenir Press, Ltd. (15s.) Toronto: 
The Ryerson Press. 

Unfortunately, the title of this book is sufficient to ensure that it wili be a commercial 
success, Public librarians will order it in hundreds without even inspecting a copy. Many 
of the public will themselves buy copies, particularly when they have read the statement 
splashed across the front cover: AT LAST THE SOVIET SCIENTISTS INVOLVED 
SPEAK OUT ON SPUTNIKS AND PROBLEMS OF THEIR CONSTRUCTION, ALSO 
ON THE NEXT GREAT STEPS TO THE MOON AND BEYOND. Of course, like most 
publisher’s blurbs, this is a mixture of half-truths and untruths. Not that we should blame 
the publishers for this—they probably know next to nothing about the subject, and in any 
case we must expect them to do everything they can to sell as many books as possible. 

A brief glance through the text or table of contents is sufficient to reveal that the major 
part of the book is not directly concerned with the sputniks, but is the usual popular account 
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of astronautics. Indeed, that is just what one would expect from the Russian title, which 
means ‘Travel into the Cosmos’’. A Russian librarian has written to say that this book 
was published by Goskul’tprosvetizdat, Moscow, in 1955, and at that date the only infor- 
mation released about the sputniks was the bald announcement that they were to be 
launched two years later. One can conclude either that the present book is a translation 
of a new edition or that some additional material has been incorporated during the trans- 
lation—and such meagre information about Sputniks I and II as does appear here contains 
nothing new. 

There is much in the English edition that can be criticized. Many of these points may 
arise from the peculiar way the text has reached us—from Russian to English, via Italian. 
For the same reason, it is impossible for the reviewer to ascribe the responsibility to any 
particular person—there are too many involved: the author, the supervisor (and, pray, why 
did Comrade Vassiliev have to be supervised ? For that matter, why did not Professor 
Dobronravov supervize him more successfully ?), the Italian translator, the English trans- 
lator, three publishers, three printers, and possibly a mixed assortment of proof-readers. It 
would obviously have been much better to have cut out as many middlemen as possible by 
going directly from Russian to English. Even if this was not feasible for some reason, and 
if it was also necessary for a translation to be made by someone unacquainted with the 
subject, it is always desirable for a translation to be read by a specialist. 

This reviewer does not feel that it is his task to list all the errors he can painstakingly 
discover (a method which is both boring to readers of the review and apt to cause unpleasant- 
ness between author and reviewer, with the editor being damned by all concerned). There 
are one or two points, however, which one cannot resist mentioning. 

In addition to the use of “‘g’’, acceleration is given in the peculiar units of feet per second. 
We also have such terms as “‘grand calories’’, “liquid combustible’, “‘body of variable 
quantity’’, ‘‘the flaming jets’’, ‘““explosions in the breach of a cannon”’ (as this dealt with a 
spacegun for lunar travel, no doubt this spelling is perfectly correct), ‘“‘thermic exchange’, 
“ratio of reverberation’’ (which appears to be the albedo of a planet), ‘“‘gas expulsion speed”’ 
(exhaust velocity). It really is time English publications stopped transliterating 
“‘Tsiolkovskii’’ into the Germanic form ‘‘Ziolkowsky’”’. 

The plates are rather a mixed bag. Some are obviously authentic Russian photographs 
(and they appear to be of too good a quality to have been copied from Vassiliev’s original 
book). Others we are apparently intended to accept as photographs, but as they are 
broadside-on views of a sputnik-launching vehicle at various stages of its flight, this seems 
improbable (the vehicle seems to be of most advanced design as it includes what I can only 
describe as a catamaran-type, multi-engined, winged booster). ‘‘A rare photograph of the 
launching of a German V. 2” is both ripe and rare—it would appear to be a small anti-tank 
rocket being fired from an armoured car. A photograph of the V.1 is included, as it is “‘one 
of the early prototypes of the rocket 'plane’’. 

One cannot recommend this edition—but it would be unfair to condemn the original 
Russian book on this evidence. I would very much like to compare the two versions, and 
in particular to see whether the author did in fact describe the V.2 as ‘‘the supreme creation 
of man’s ingenuity’. That seems a peculiar statement to come from one of the “Soviet 
scientists involved’”’ in sputnik construction—one would not even expect it from Dr. 
von Braun and his co-workers at Peenemiinde, who would have had much more excuse 
for making it. 

G. V. E. THoMPsoN. 


Digital Differential Analysis: An Applications Manual for Digital and Bush Type 
Differential Analyzers. By George F. Forbes. Fourth Edition. 10} x 84 in. 
Pp. (xiv) + 154 + (xxxviii), illustrated. 1957. Pacoima, California: George F. Forbes, 
10117, Bartee Avenue. ($5.00.) 

This is a manual for mathematicians, engineers and researchers who intend to use a 
digital differential analyzer (DDA) or a Bush analyzer. It briefly describes the mechanisms 
of the DDA, which, as its name implies, is a digital (pulse) computer employing vacuum 
tubes, crystal diodes and magnetic drum memory. The main part of the manual is devoted 
to a discussion of programming problems for the machine, and it covers a wide variety of 
mathematical functions and techniques. 

Of particular interest to readers of this Journal is Chapter XII, which discusses the 
trajectory of a powered rocket in the presence of a central force. Papers by Forbes on this 
subject have appeared in this Journal (1950, 9, 75; 1951, 10, 194) and a summary of Chapter 
XII appeared in ].B.J.S., 1955, 14,85. Detailed solutions of this problem are not givenin the 
present work, since it is the method of solution using the computer that is considered 
important for this manual, and rightly so. 

D. J. CASHMORE. 
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Rockets, Missiles, and Space Travel. By Willy Ley. Revised edition. 8$ 

Pp. xv + 528, with 85 figures and 31 plates. 1957. New York: Viking Press ( $6.75.) 

London: Chapman and Hall, Ltd. (£2 10s.) 

Good wine needs no bush and this work needs no recommendation from me. It has 
long held a place as the bible of rocketry, presenting as it does an authoritative history of 
the subject, together with some prophetic chapters about the future! Originally appearing 
as Rockets in May, 1944, it has steadily grown in size to keep pace with the rapid post-war 
developments in astronautics. After three printings, it reappeared in February, 1947, 
under the title Rockets and Space Travel (two printings). Then in June, 1951, a new 
edition with the present title was issued; this went through six printings. 

While chapters dealing with events prior to 1930 are substantially unchanged, this 
new edition has been extensively revised and enlarged (since the 1951 edition, the size has 
increased by 92 pages) to cover the American upper atmosphere research and Vanguard 
programmes, new missiles, etc. The story of Peenemiinde has been completely rewritten. 
Many additions have been made to the valuable appendices and bibliography. 

Mr. Ley has long held a place as the premier historian of astronautics and possession 
of this book is an essential for anyone with a serious interest in the subject. Of course, in 
common with all other astronautical books it is already out of date (there is naturally no 
mention of the Russian sputniks), but even a loose-leaf service would find it hard to keep 
pace with developments in this dynamic, rapidly-developing branch of science. However, 
no doubt, Mr. Ley is already hard at work at further revisions and there is much to be said 
for the acquisition of each new edition as it appears—they are themselves an indication of the 
progress that is being made. In the fourteen years that have elapsed since Rockets made 
its appearance many of the future possibilities it mentioned have become realities. 

G. V. E. THompson. 
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